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RESEARCH SUMMARY 


Data on regeneration of 10 conifer species commonly 
found in the Northern Rocky Mountains of the Western 
United States were collected by siratified random sam- 
pling of 12,128 '/s00-acre plots in 537 stands. These 
data were collected in Montana, central Idaho, and north- 
ern Idaho, representing five habitat type series having 
various combinations of regeneration harvest methods 
and site preparations. Also included is the influence of 
western spruce budworm (Choristoneura occidentalis) 
on regeneration success. The resulting Regeneration 
Establishment Model is a part of the Prognosis Model. 
Together, these models can be used to predict growth 
and yield for even-age and uneven-age management 
systems. The regeneration model predicts the prob- 
ability of stocking, seedling density, species composi- 
tion, and seedling heights 2 to 20 years after harvest 
for various combinations of regeneration methods, site 
preparations, habitat types, and budworm defoliation 
histories. 
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Predicting Regeneration 
Establishment Witn the 
Prognosis Model 


Dennis E. Ferguson 
Clinton E. Carlson 


INTRODUCTION 


Forests of Idaho and Montana are quite diverse. Environmental conditions 
support forest types ranging from dry ponderosa pine (Pinus ponderosa) for- 
ests, to near-rain forests dominated by western redcedar (Thuja plicata) and 
western hemlock (Tsuga heterophylla), to high-elevation subalpine fir (Abies 
lasiocarpa) and whitebark pine (Pinus albicaulis) forests where the growing 
season is short. These forests are highly valued for timber, wildlife, water, 
grazing, recreation, biological diversity, and esthetics. 

Management often involves the manipulation of conifers through timber 
harvests, natural and artificial regeneration, site preparation, weedings, 
and thinnings; control of competing vegetation such as shrubs, forbs, and 
grasses; and control of insects, diseases, and wildfires. Manipulation of for- 
ests changes the course of secondary succession which, in turn, affects the 
future of these ecosystems. 

Land managers need to know how their actions affect the course of second- 
ary succession. One particularly important effect is the regeneration of coni- 
fers following timber harvest. Can cutover areas be regenerated in a timely 
manner to desirable species having acceptable growth rates? Which of sev- 
eral alternative silvicultural prescriptions should be used to reach desired 
future conditions? 

Land managers compare alternative silvicultural prescriptions to determine 
which one best meets management objectives. These comparisons need to 
be quantified in a biologically meaningful and statistically sound manner. 
Fortunately, the necessary information is available to quantify the regener- 
ation process in much of the Western United States. 

The first requirement for predicting regeneration success is an understand- 
ing of the silvical characteristics of conifer species. Many research reports deal 
with silvics of conifers found in the interior Western United States. For exam- 
ple, Haig and others (1941) summarize many years of research in northern 
Idaho that provides information on silvical characteristics of conifers and 
the effects of silvicultural treatments on regeneration success. Many other 
studies are also available; a good place to begin is Burns and Honkala (1990). 

The second requirement for predicting regeneration success is an ecologi- 
cal land classification system. The habitat type system, first developed by 
Daubenmire (1952), provides the framework for modeling secondary succession. 
Today, more detailed habitat type classifications are available for Montana 
(Pfister and others 1977), central Idaho (Steele and others 1981), and north- 
ern Idaho (Cooper and others 1991). 

The third requirement is a large number of managed stands that can be sam- 
pled to represent independent variables in a regeneration model—number 
of years since disturbance, site preparation, residual overstory density and 


1 


species composition, habitat type, and so on. The Northern Region (Region 1) 
and Intermountain Region (Region 4) of the Forest Service, U.S. Department of 
Agriculture, have a large number of managed stands that can be retrospectively 
examined because there are adequate historical records. Additional stands 
with adequate historical records are available from other major ownerships, 
such as Potlatch Corporation, Idaho Department of Lands, Montana Depart- 
ment of Forestry, Champion Corporation, or the Bureau of Indian Affairs 
and Bureau of Land Management in the U.S. Department of the Interior. 

These three factors enabled a retrospective examination of regeneration 
success in cutover stands and the development of a quantitative regenera- 
tion model. The Regeneration Establishment Model described in this paper 
is a part of the Prognosis Model (Stage 1973; Wykoff and others 1982). The 
Regeneration Establishment Model predicts establishment of new trees that 
are added to the inventory of trees being projected by the Prognosis Model. 
The Prognosis Model implements user-specified management prescriptions 
and predicts accretion and mortality for the stand. The Prognosis Model with 
the Regeneration Establishment Model can be useful to managers in the fol- 
lowing situations: 


1. To decide which of several prescriptions best meets reforestation objectives. 

2. To link regeneration to predictions of long-term growth and yield. 

3. To predict regeneration in long-term projections of succession or uneven- 
age management. 

4. To help develop harvesting schedules. 

5. To help in other areas of forest planning where regeneration is important 
for watershed, wildlife, esthetics, and so on. 


Version 1 of the Regeneration Establishment Model was developed for the 
grand fir-cedar-hemlock ecosystem in northern Idaho and adjacent portions 
of Montana and Washington (Ferguson and others 1986). Version 2 of the 
model, reported in this paper, expands the geographic area, includes additional 
habitat types, and includes the influence of western spruce budworm (Chor- 
istoneura occidentalis) on regeneration success. 

Version 2 of the Regeneration Establishment Model is calibrated for west- 
ern Montana as well as central and northern Idaho. Figure 1 shows the loca- 
tions of sampled stands. Habitat type series included in the model are the 
Pseudotsuga menziesii, Abies grandis, Thuja plicata, Tsuga heterophylla, 
Abies lasiocarpa, and some of the Tsuga mertensiana series. 

The effects of western spruce budworm defoliation on regeneration success 
are represented in the regeneration model. Western spruce budworm reduces 
growth, crown volume, and cone crops of host trees, and increases top-kill, 
mortality, deformities, and predisposition to other insects and disease. Host 
species in the study area are western larch (Larix occidentalis), Douglas-fir 
(Pseudotsuga menziesii), grand fir (Abies grandis), Engelmann spruce (Picea 
engelmannii), and subalpine fir (Abies lasiocarpa). 

To be consistent, we use scientific names and four-letter abbreviations 
when discussing habitat types. An example is the Abies grandis/Clintonia 
uniflora habitat type, abbreviated ABGR/CLUN. When discussing tree spe- 
cies, we use the common name; for instance, grand fir. 

This paper presents details about data collection, predictive equations, model 
development, and model behavior. A user’s guide to version 2 of the Regen- 
eration Establishment Model is available (Ferguson and Crookston 1991). 
Silvicultural interpretations of the model will be presented in another paper 
(Carlson and Ferguson in preparation). 
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Figure 1—Location of study sites for version 2 of the Regeneration Establishment 
Model. Triangles are stands sampled by Carlson in 1979-82. Squares are stands 
sampled by Ferguson and others (1986) in 1975-76. Circles are stands sampled 
by Ferguson in 1979-82. 


STUDY DESIGN 


Data for version 2 of the regeneration model were collected during three 
studies. The grand fir-cedar-hemlock ecosystem of northern Idaho, sampled 
in 1975-76, was used to develop version 1 of the regeneration model (Ferguson 
and others 1986). Data from western Montana were collected by Carlson from 
1979-82. Data in central Idaho and the Pseudotsuga menziesii and Abies 
lasiocarpa series in northern Idaho were collected by Ferguson from 1979-82. 
Data collection was similar among the three studies. These data were com- 
bined to develop version 2 of the regeneration model. 

A stratified random design was used to select stands that were harvested 
2 to 20 years prior to sampling. Stratification variables were habitat type 
series, site preparation, regeneration method, and history of budworm defolia- 
tion. Four to five stands within each stratification category were randomly 
chosen for sampling. 


Although stands were used to select study areas, we were interested in 
sampling microsites within stands. A series of Y400-acre plots were sampled 
in each stand. After obtaining stand histories and stereo aerial photographs 
for each selected stand, either transects or a sampling grid was used to sys- 
tematically locate about 25 plot centers per stand. In Idaho, transects crossed 
through areas observable on aerial photographs that varied in such things 
as aspect, overstory density, site preparation, and topographic position. In 
Montana, a sampling grid was used to locate plot centers. 

Each sample point was the center of a Ysoo-acre circular fixed area plot and 
a variable radius plot for sampling the overstory (if present) with a 10 basal 
area factor prism. Each plot was made as independent as possible by recording 
plot attributes—habitat type, slope, aspect, type of site preparation, topographic 
position, and residual overstory density by species. Systematically spacing 
plots also helped maintain independence among plots within stands because 
plots were dispersed throughout the stand. The year of cutting, geographic 
location, elevation, and budworm defoliation history were stand variables 
that did not differ among plots within a stand. 

Plots within stands could represent different “years since last disturbance” 
due to different years of site preparation. Untreated plots were last disturbed 
when the stand was harvested. Plots mechanically prepared or burned were 
last disturbed when the site preparation took place—either, or both, could 
occur after the harvest year. Because the year of site preparation can take 
place after the year of harvest, we use “years since last disturbance,” rather 
than stand age, to predict regeneration. 

Established conifer regeneration on each plot in the Idaho studies was 
counted by species; in the Montana study up to 10 established trees were re- 
corded per plot. Species included in the model are shown in table 1. Mountain 
hemlock (T'suga mertensiana) could not be included in the model because of 
insufficient data. Table 1 also lists the minimum establishment heights by 
species—these are 0.5 foot for shade-tolerant species and 1.0 foot for shade- 
intolerant species. Minimum heights roughly correspond to the height of a 
3-year-old tree. Maximum size for all regeneration was 2.95 inches diameter at 
breast height (d.b.h.). Tree heights were measured as of the end of the pre- 
vious growing season to ensure that all trees sampled during one field season 
used the same standards for establishment. 

Tree height to the nearest 0.1 foot, tree age at groundline, and tree condi- 
tion (damage, insects, disease, and so on) were recorded on a subsample of 
regeneration called “best trees.” Best trees were chosen by the following rules: 


1. Select the two tallest trees on each Yso0-acre plot regardless of species. 

2. Select the one tallest tree of each additional species present on the plot. 

3. If the first two rules do not total four trees, select in order of descending 
height from any remaining trees until at least four trees, if present, are chosen. 


One site preparation treatment was recorded for each Ysoo-acre plot. It was 
the treatment that covered the largest proportion of the plot—mechanical, 
burn, no treatment, or road. The “road” category included all road cuts, road 
fills, and unmaintained roadbeds. Plots that occurred on maintained road- 
beds were not sampled. The intensity of each site preparation was not re- 
corded because this was more difficult to determine as years since disturbance 
increased. 

National Forests were used as independent variables in data analyses be- 
cause they divide the study area into units that can account for variation in 
latitude, longitude, climate, and so on. 


Table 1—Scientific name, common name, abbreviation, and minimum establishment height for 
regeneration. Species are listed in approximate order of increasing shade tolerance. 
This order will be maintained in tables that provide equations by species 


Prognosis 
Model Minimum 
Scientific name Common name abbreviation height 
Feet 

Pinus ponderosa Doudgl. ponderosa pine PP 1.0 
ex Laws. var. 
ponderosa 

Larix occidentalis western larch L 1.0 
Nutt. 

Pinus contorta Dougl. lodgepole pine LP 1.0 
ex Loud. 

Picea engelmannii Engelmann spruce S ifs) 
Parry ex Engelm. 

Pseudotsuga menziesii Douglas-fir DF 1.0 
var. glauca (Beissn.) 
Franco 

Pinus monticola Dougl. western white pine WP 1.0 
ex D. Don 

Abies grandis (Doug. grand fir GF As} 
ex D. Don) Lindl. 

Tsuga heterophylla western hemlock WH 5 
(Raf.) Sarg. 

Abies lasiocarpa subalpine fir AF Ao) 
(Hook.) Nutt. 

Thuja plicata Donn ex western redcedar Cc 5 


D. Don 


The presence of shrubs, forbs, and grasses was recorded on each ¥300-acre 
plot. Each species was characterized by recording average height and per- 
centage of the plot area covered by a vertical projection of the crown on the 
ground. Coverage could exceed 100 percent because of overlapping layers 
of shrubs, forbs, and grasses. These data were not used to predict regenera- 
tion success because it was not possible to retrospectively determine cover- 
age at the time of seedling establishment. A model of secondary succession 
for shrubs, forbs, and grasses has been developed using the Idaho data 


(Moeur 1985). 


MODELING TECHNIQUE 


The technique for developing equations follows the example of Hamilton and 
Brickell (1983) for two-state systems. Each plot in the sample is in one of 
two states: stocked with at least one established seedling or nonstocked. All 
plots are used to develop equations predicting the probability of a plot being 
stocked. Then, only stocked plots are used to estimate the number of trees 
on the plot, number of species, species composition, and seedling heights. 

The probability of stocking ¢ years after disturbance is estimated by a logistic 
equation (Hamilton 1974) within the interval [0,1]. After the attributes of 


Logistic 
Regression 
Equations 


stocked plots are estimated, the probability of stocking is used to scale stocked 
plot attributes to a per-acre basis. As an example, suppose five plots are being 
projected and the following are predicted for the probability of stocking (PS) 
and number of trees per stocked plot (TPSP): 


Predicted Predicted number Predicted trees 
probability of trees per per acre = 
of stocking stocked plot (PS * TPSP * 300)/N 
(PS) (TPSP) where N is 5 plots 
0.60 3 108 
25 1 15 
.80 2 96 
.50 5 150 
45 Ht 2 
Predicted total 
trees per acre 396 


The third column shows how the probability of stocking is used to scale results 
of stocked plot analyses to a per-acre basis. The “300” in the third column 
is the inverse of the plot size used in the study and the “N” is the number of 
plots being used to represent the stand. Predicted total trees per acre is 396. 
Species, size, and other tree attributes are determined during other steps 
in the model. 


Many of the regeneration data are distributed dichotomously—plots are 
either stocked with at least one established seedling or they are nonstocked, 
each of the 10 species was either established on a plot or it was not, and so 
on. Response surfaces are sigmoid-shaped. 

A nonlinear logistic algorithm called RISK (Hamilton 1974) was used to 
predict probabilities for dichotomously distributed dependent variables. 
The logistic equation is: 


P = (1+e6%)) (1) 


where 
P = probability 
e = the base of natural logarithms 
B; = vector of regression coefficients 
X; = vector of independent variables. 


The predicted probability (P) is continuous and bounded within the interval 
[0,1]. RISK can use up to 30 independent variables. 

Goodness of fit for equations was evaluated at the 0.05 significance level. 
RISK output includes four types of statistics that can be used to evaluate 
goodness of fit: Student-t ratios for each regression coefficient, chi-square 
values for 21 divisions within the probability interval [0,1], total chi-square 
for the 21 divisions, and an analysis of variance table showing the amount 
of variation explained by regression. 

A fifth goodness-of-fit criterion was how close the error mean square was to 
1.0 in the analysis of variance table. As the population size increases, the limit- 
ing value for the error mean square is 1.0, rather than 0.0 (Hamilton 1974). 
We found that in successive iterations of the RISK algorithm, the error mean 
square would approach 1.0 unless the equation was unstable; thus, error 
mean squares substantially different from 1.0 usually indicate poor equa- 
tion formulation. 


Linear Regression 
Equations 


Weibull 
Distributions 


Habitat Type 
Groupings 


Linear regression was used to estimate heights of best trees. Goodness of 
fit was evaluated at the 0.05 significance level. The best transformation of 
the dependent variable was found to be the natural logarithm of tree height. 
Various transformations of independent variables were also explored and used. 


In some steps of the regeneration model, it is important to mimic the dis- 
tribution of the dependent variable. An example is the number of seedlings 
established on stocked plots, which has a reversed J-shape distribution. For 
these situations, the distribution is modeled with Weibull cumulative den- 
sity distributions (Bailey and Dell 1973). The form of the equation for a two 
parameter Weibull cumulative density distribution is: 


F(x) = 1 -EXP[-(x/B)] (2) 
where 


F(x) = cumulative density 
x = dependent variable being modeled 
B = Weibull scale parameter 
C = Weibull shape parameter. 


The area bounded by the Weibull equation is within the interval [0,1]. 
Equation 2 can be solved for x as shown by Bailey and Dell (1973): 


x = B[-In(1-F(x))]“" (3) 


A uniformly distributed pseudo-random number in the interval [0,1] can 
be substituted for F(x) in equation 3, resulting in an unbiased choice for x. 
Plot replication in the regeneration model results in enough random draws 
that the distribution is approximated. 

Bailey and Dell (1973) also discuss coefficients for the C (shape) parameter. 
If C is less than 1, the distribution has a reversed J shape. If C equals 1, 
the distribution is exponential. If C is between 1 and 3.6, the distribution is 
mound-shaped and positively skewed. When C is approximately 3.6, a nor- 
mal distribution results. And when C is greater than 3.6, the distribution is 
mound-shaped and negatively skewed. 


Version 1 of the regeneration model predicted regeneration on four 
broadly defined habitat types using the habitat type classification system of 
Daubenmire and Daubenmire (1968). Since then, classification of habitat 
types has been refined for northern Idaho (Cooper and others 1991). We were 
fortunate that Neiman, one of the authors of the 1991 northern Idaho classi- 
fication, was available to update data collected during 1975-76 to the 1991 
classification. Each plot sampled during the 1975-76 study was reclassified 
by reviewing the original field sheets on which field crews had listed species 
of shrubs, forbs, and grasses by plot. 

Thus, for northern Idaho we used the habitat type classification of Cooper 
and others (1991), for Montana we used the classification of Pfister and others 
(1977), and for central Idaho we used the classification of Steele and others 
(1981). This resulted in 93 different habitat types and phases being sampled. 

It was necessary to combine similar habitat types into groups to reduce 
independent variables to a reasonable number. This was done before fitting 
equations. Groupings were made after statistical analyses by habitat types 
that compared the percentage of stocked plots, trees per stocked plot, and 
species per stocked plot. These analyses provided ideas for preliminary 
groupings. Preliminary groupings were reviewed by forest ecologists, and 
their suggestions were used to develop the final groupings. 
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Western Spruce 
Budworm 
Defoliation 
History 


Planting 


Variables in the 
Model 


The resulting 16 habitat type groups are shown in table 2. Analyses used 
to develop the groups are found in appendix A. The groups combine habitat 
types within climax series that are similar during the regeneration phase of 
stand development. Habitat type groups would likely be different for other 
purposes—diameter growth, mortality, canopy cover, and so on. 


The influence of western spruce budworm on regeneration success was a 
major focus of the studies conducted in Montana and Idaho during 1979-82. 
Defoliation histories were collected by stand for the period 5 years before 
the harvest to the year prior to field sampling. 

In Montana, past budworm activity was indexed by analysis of increment 
cores using the procedure of Carlson and McCaughey (1982). In Idaho, his- 
torical aerial defoliation maps were used to determine which years, if any, 
a stand was within a defoliated area. Using aerial defoliation maps, it was 
also possible to determine defoliation histories for stands in the 1975-76 
study. The vast majority of stands from the 1975-76 study were not in areas 
defoliated by budworm. 

The historical aerial defoliation maps outline areas of low, moderate, and 
heavy defoliation as seen by an entomologist from an airplane. Boundaries 
between levels of budworm defoliation are approximate. There is also varia- 
tion in defoliation within outlined levels of defoliation. We felt the appropri- 
ate way to use budworm histories in developing a regeneration model was 
as a dichotomous variable. For each year, the stand either was in a defoli- 
ated area or it was not. This approach also allowed consistent use of data 
collected in Montana and Idaho. 


Both planted and unplanted stands were sampled. Sampling only naturally 
regenerated stands may have biased the sample, because stands that were 
difficult to regenerate were more likely to have been planted. During the 
first few years of data collection, crews attempted to identify planted trees on 
plots, but this was not easy to determine retrospectively. Plots from planted 
stands were coded for analyses as a class (dummy) variable. Coefficients 
estimated by regression analysis for planted stands are not used in the model, 
so that predictions are only for natural regeneration. Model users then specify 
planting for the stand, which is added to the natural regeneration predicted 
by the model. They select the species, density, year of planting, and survival 
of planted trees. 


The variables used in the regeneration model are those generally recorded in 
forest inventories in the Northern Rocky Mountains. This means that most 
inventories are already designed to collect the information needed to use the 
regeneration model, or could easily do so. 

Table 3 provides definitions for variables used in the regeneration model. 
These names for variables are used in tables of regression coefficients 
(appendix B). 

During equation development, goodness of fit was evaluated at the 0.05 
significance level. Nonsignificant independent variables were eliminated 
from equations except those we considered to be part of groups. An example 
of a group is class variables for site preparations. The entire group was 
kept in the equation if one of the coefficients was significant. 


Table 2—Number of plots by habitat type and habitat type group. 
Within a habitat type group, habitat types are listed in de- 
creasing order of number of plots sampled. Subsequent 
tables will use only the first habitat type to identify the group 


Habitat types’ No. plots 
Group (number of '/300-acre plots) in group 
1 PSME/VAGL, LIBO, VACA 296 
(178) (92) (26) 
2 PSME/CARU, CAGE, AGSP, FEID 417 
(264) (122) (26) (5) 
3 PSME/PHMA, ACGL 832 
(815) (17) 
4 PSME/SYAL, SPBE, SYOR, ARUV, Misc. 887 
(568) (223) (75);* (18) (3) 
5 ABGR/LIBO, CLUN-XETE 273 
(199) (74) 
6 ABGR/XETE, VAGL, COOC, VACA 451 
(248) (i777) een (3) an (3) 
ABGR/CLUN (except CLUN-XETE) 1,865 
ABGR/SPBE, ACGL, PHMA, ASCA, SETR 989 
(447) (251) (GS) “(38)) (2) 
9 THPL/All 2,180 
10 TSHE/AIl 1,387 
11 ABLA/VAGL, VASC, VACA 342 
(176) (110) (56) 

12 ABLA/XETE, LIBO 675 
(533) (142) 

13 ABLA/CLUN, GATR 680 
(673) (7) 

14 ABLA/CAGE, CARU, ACGL, SPBE 228 
(105) (73) (36) (14) 

15 ABLA/MEFE, ALSI 
(353) (31) 

TSME/CLUN, XETE, MEFE, STAM 460 

(31) (Sil) ee (9) ea (9) 

16 ABLA/CACA, STAM, LUHI 166 
(107) (55) (4) 

Total 12,128 


‘Species codes for habitat type abbreviations: 


ABGR_ Abies grandis ABLA Abies lasiocarpa 

ACGL Acer glabrum AGSP_ Agropyron spicatum 

ALS!| Alnus sinuata ARUV Arctostaphylos uva-ursi 
ASCA Asarum caudatum CACA Calamagrostis canadensis 
CAGE Carex geyeri CARU Calamagrostis rubescens 
CLUN  Clintonia uniflora COOC Coptis occidentalis 

FEID Festuca idahoensis GATR_ Galium triflorum 


LIBO __Linnaea borealis LUHI 


MEFE Menziesia ferruginea 
PSME Pseudotsuga menziesii 
SPBE Spiraea betulifolia 
SYAL Symphoricarpos albus 


PHMA 
SETR 
STAM 
SYOR 


Luzula hitchcockii 
Physocarpus malvaceus 
Senecio triangularis 
Streptopus amplexifolius 
Symphoricarpos oreophilus 


THPL Thuja plicata TSHE Tsuga heterophylla 
TSME Tsuga mertensiana VACA _ Vaccinium caespitosum 
VAGL_ Vaccinium globulare VASC_ Vaccinium scoparium 


XETE Xerophyllum tenax 


Table 3—Definition of variables used in the Regeneration Establishment Model 


Variable 
ASP 
SLO 
TIME 
REGT 
BWAF 


BWB4 
SQREGT 


SQBWAF 


ELEV 
BAA 
NONE 
MECH 
BURN 
ROAD 


PLANT 


TPSP 


OVER 
BOTTOM 
LOWER 
MIDSLOPE 
UPPER 
RIDGE 
AGE 


Definition 
Plot aspect converted to radians. 
Plot slope tangent (slope percent/100). 
Number of years since last disturbance to plot. 
Number of years since last disturbance without budworm. 


Number of years since last disturbance with budworm. 
Note: REGT + BWAF = TIME 
Number of budworm years in the 5 years before last disturbance. 


Sum of square roots for each year without budworm. 
SQREGT = X{(TIME") — [(TIME-—1)'7]} for those years without budworm. 
Suppose there was no budworm in years 1, 2, and 5. 


Then: 
YEAR = 1 2 3 4 5 
TIME*?— [(TIME—1)"”)] = 1.000 0.414 0.236 


SQREGT = 1.000 + 0.414 + 0.236 = 1.650 


Sum of square roots for each year with budworm. 
SQBWAF = 2{(TIME") — [(TIME—1)"?]} for those years with budworm. 
Suppose there was budworm in years 3 and 4. Then: 


YEAR = 1 2 3 4 5 
TIME*?— [(TIME-1)"”)] = 0.318 0.268 
SQBWAF = 0.318 + 0.268 = 0.586 


The effect of this variable is that budworm defoliation near the beginning 
of the regeneration period will be more detrimental than defoliation near 
the end of the regeneration period. 

Note: (SQREGT + SQBWAF)? = TIME 


Stand elevation to the nearest 100 feet; for example, 35 = 3,500 feet. 
Plot residual overstory basal area (ft?/acre), reset to 1.0 if 0.0. 

Class variable for no site preparation to the plot. 

Class variable for mechanical disturbance to the plot. 

Class variable for burn disturbance to the plot. 


Class variable for plots that occur on road cuts, road fills, or unmaintained 
roadbeds. 


Class variable for planted stands (probability of stocking equations and 
equations for determining the number of species on stocked plots) or the 
stand planted to the same species being predicted (probability of species 
occurrence equations). 


Number of regeneration-size trees on the plot, conditional on the plot 
being stocked with at least one established seedling. 


Class variable for the presence of the same species in the variable radius plot. 
Class variable for bottom topographic position. 

Class variable for lower-slope topographic position. 

Class variable for midslope topographic position. 

Class variable for upper-slope topographic position. 

Class variable for ridge topographic position. 

Tree age at groundline. 
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RESULTS 


This section describes equations that were statistically derived from the 
data. Table 4 shows the range of data for independent variables used to 
develop the model. Data represent 2-20 years since disturbance, all aspects, 
0-390 square feet of residual overstory basal area, 2,400-7,400 feet elevation, 
four site preparation methods, 0-110 percent slopes, and up to 16 years of 
budworm defoliation. The year of harvest for stands included in this study 
ranged from 1959 through 1978. 

Table 5 shows the number of plots sampled by classification variables of 
regeneration method, site preparation, habitat type series, and budworm 
history. Except for the road site preparation category, this is the same in- 
formation matrix used in the stratified random selection of stands. Some 


Table 4—Number of plots by classes for independent variables used to develop version 2 of the Regeneration Estab- 
lishment Model. Classes used below show the range of data—analyses used the data as continuous vari- 
ables unless otherwise noted in appendix B 


Class No. plots Class No. plots Class No. plots 
Years since Site preparation Aspect (degrees) 
last disturbance None 5,648 Flat 204 
2 58 Mechanical 3,055 346-45 1,763 
3 768 Burn 2,742 46-105 2,009 
4 849 Road 683 106-165 2,171 
5 974 Regeneration 166-225 1,996 
6 986 method 226-285 2,016 
7 934 286-345 1,969 
8 922 Clearcut 4,884 ; 
9 914 Seedtree 2,207 Topographic 
10 762 Shelterwood 2,872 position 
Selection 2,165 Bottom slope 628 
Note The) Lower slope 820 
13,14 1,233 Regeneration mode Midslope 7.720 
15,16 1,062 Natural 7,624 Wanercione 2 634 
17,18 631 Planted 4,504 suse 
419.20 586 Ridgetop 326 
: Slope (percent) No. of budworm 
Residual basal Flat 313 ines he 5 
area (ft?/acre) 1-10 2,187 heron he . is eas 
0 6,104 11-20 2,822 Senet eo 
10 1,089 21-30 2,637 1 68 4 
20 1,034 31-40 1,794 2 579 
30-40 1,344 41-50 1,220 3 267 
50-60 913 51-80 1,095 4 384 
70-80 609 81-110 60 
5 860 
peg’ ce Nearest National 
110-200 549 Forest No. of budworm 
210-390 45 : years after harvest 
Gallatin 178 0 7232 
Elevation Helena 371 1 338 
(nearest 100 feet) Deerlodge 321 2 588 
24-27 204 Bitterroot 461 3 678 
28-31 852 Flathead 199 4 635 
32-35 997 Lolo 1,446 5 474 
36-39 1,030 Kootenai 356 6 388 
40-43 1,879 Panhandle 2,933 7 261 
44-47 1,668 Colville 271 8 401 
48-51 1,721 Clearwater 2,527 9 366 
52-55 1,178 Nez Perce 1,037 10 318 
56-59 1,079 Payette 1,399 14 132 
60-63 741 Boise 629 12 171 
eeu can 13,14 100 
68-74 334 , 


15,16 46 
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Table S—Number of '00-acre plots sampled by classification variables 


Budworm __ Habitat Regeneration method and site preparation 


defoliation type Clearcut Seedtree Shelterwood Selection 
history series' None Mech Burn Road None Mech Burn Road None Mech Burn Road None Mech Burn Road 
No PSME 80) awl 705 ail 82 30n eet 7 176 124 61 18 302 67 2 13 


ABGR 197 Sea 63 74 229. 1389 87 41 135 ~=«61 19 3 211 53 18 9 
THPL 193. 193 251 44 133 129 4 31 265 31 10 29 256 31 3 30 
TSHE 275m, W230 232hg 53 150 34 37 16 V2 OS 2 13 122 48 35 17 


ABLA 149 121 218 18 134 109 51 14 139 68 20 13 106 16 5 8 
Subtotals 894 679 934 167 728 441 156 109 887 335 112 76 997 215 86 77 
Yes PSME 66 155 146 8 138 63 75 3 275 ~—s «61 72 3 153 50 8616 0 
ABGR 174 251 493 51 99 560 4700 pil 304 171 HO 32 113 90) 22) 13 
THPL 1220 24p Set Ot 16 29 16 0 2 99° 35) «14 18 53 10 2 6 
TSHE 2 1 U 1 0 0 0 0 8 6 5 2 8 1 0 0 
ABLA 94 168 282 48 89 91 40 10 135 89 49 14 181 47 13 12 
Subtotals 458 599 1,029 124 355 226 162 30 821) 362) 7210)) 169 508 198 53 31 
Totals 1,352 1,278 1,963 291 1,083 667 318 139 1,708 697 322 145 1,505 413 139 108 

Y= 12,128 


'See table 2, footnote 1, for habitat type abbreviations. 


combinations were difficult to find; for example, selection cuttings that were 
burned. Also, spruce budworm defoliation was rarely found in Tsuga hetero- 


phylla habitat types. 
Probability of Probability of stocking equations predict the probability of at least one 
Stocking established seedling on a ¥00-acre plot. Separate equations were developed for 


the Pseudotsuga menziesii, Abies grandis, Thuja plicata/Tsuga heterophylla, 
and Abies lasiocarpa series (appendix B, table 9). A fifth equation was devel- 
oped for the 683 plots that were sampled on road cuts, road fills, and unmain- 
tained roadbeds. 

The equations given in appendix B, table 9, were exercised to show the 
effect of important independent variables on the probability of stocking. In 
figures 2 to 9, independent variables are held constant except those being 
displayed. Default conditions are an Abies grandis/Clintonia uniflora habi- 
tat type on the Idaho Panhandle National Forests, north aspect, 30 percent 
slope, 10 years since last disturbance, 3,500 feet elevation, no residual over- 
story, no site preparation, and no planting. 

Figure 2 illustrates the interaction of slope and aspect on predicted prob- 
ability of stocking. On southerly slopes, the probability of stocking decreases 
quite rapidly with increasing slopes. On northerly aspects, the probability 
of stocking increases only slightly as slope increases. 

Increases in the probability of stocking over time differ by aspect (fig. 3). 
North aspects increase the most rapidly and attain the highest stocking level. 
South aspects have the poorest stocking over time. East and west aspects 
are intermediate between north and south. For equations predicting the 
probability of stocking, the highest rate of stocking over time occurs at the 
optimum aspect and the lowest rate of stocking occurs at the poorest aspect. 

Habitat types also differ in the rate of stocking over time. These differences 
are illustrated in figure 4. Of the four habitat types shown (each representing 
a different equation), the Thuja/Tsuga series has the highest rate of stocking 
and achieves the highest stocking level. The probability of stocking on the 
Abies grandis series increases nearly as quickly as the Thuja/Tsuga series, 
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Figure 2—The effect of slope and aspect on predicted probability 
of stocking for an Abies grandis/Clintonia uniflora habitat type. 


Probability of Stocking 


315 


0 5 10 15 
Years Since Last Disturbance 


Figure 3—The effect of aspect and time on predicted probability 
of stocking for an Abies grandis/Clintonia uniflora habitat type. 
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Figure 4—The effect of habitat type and time on predicted probability 
of stocking. 


but never achieves as high a stocking level. The Abies lasiocarpa series has 
a lower rate of stocking than the Thuja/Tsuga and Abies grandis series and 
also does not achieve as high a stocking level. The Pseudotsuga menziesii 
series has the lowest stocking level and the lowest increases over time of all 
the series. 

Elevation is another important independent variable. For the Abies grandis 
and Thuja/Tsuga series, the effect of elevation on the probability of stocking 
is quadratic (fig. 5). In the Thuja/Tsuga series, the optimum elevation is 
about 3,500 feet. The optimum elevation in the Abies grandis series is near 
4,500 feet. The quadratic effect of elevation in the Abies grandis series is less 
than that in the Thuja/Tsuga series. In the Pseudotsuga menziesii series, 
increasing elevation decreases the probability of stocking. Elevation was 
not an important independent variable in the Abies lasiocarpa series. 

The effect of residual overstory basal area on the probability of stocking 
changes with aspect and slope. The general trend of residual basal area on 
the probability of stocking for the Abies grandis and Thuja/Tsuga series is 
quadratic. On north aspects, optimum basal area decreases as slope increases 
(fig. 6). This means that-a steep north-facing slope does not need as much 
shelter as a gentle north-facing slope. On south aspects, optimum basal area 
increases as slope increases (fig. 7). Steep south-facing slopes need more 
shelter than gentle south-facing slopes. 

Residual basal area was not important in the equation for the Pseudotsuga 
menziesii series. In the Abies lasiocarpa series, the natural log of residual 
basal area was significant, meaning that the probability of stocking increases 
rapidly as basal area increases at low densities, but there is little additional 
effect as basal area gets progressively more dense. 
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Figure 5—The effect of elevation and habitat type on predicted 
probability of stocking. 
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Figure 6—The effect of residual overstory basal area and slope 
for a north aspect on predicted probability of stocking for an Abies 
grandis/Clintonia uniflora habitat type. 
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Figure 7—The effect of residual overstory basal area and slope 
for a south aspect on predicted probability of stocking for an 
Abies grandis/Clintonia uniflora habitat type. 


The effect of western spruce budworm on the probability of stocking is not 
very dramatic. This is not surprising since nonhost species stock plots as do 
budworm host species. Lower rates of stocking over time occurred in the 
Pseudotsuga menziesii, Thuja/Tsuga, and Abies lasiocarpa series, but not in 
the Abies grandis series. Additionally, budworm defoliation in the 5 years 
prior to harvest significantly decreased the probability of stocking in the 
Abies lasiocarpa series. Figure 8 illustrates the probability of stocking over 
time on the Pseudotsuga menziesii series. The upper line has no budworm. 
The lower line illustrates the probability of stocking if budworm defoliation 
occurred every year. After 20 years, the difference in the probability of 
stocking is 0.07. 

Coefficients for spruce budworm effects in appendix B are sometimes not 
significant at the 0.05 level. Budworm defoliation reduces the rate of stock- 
ing of host species; it does not prevent stocking. For example, in the equation 
for the probability of stocking for Pseudotsuga menziesii series (appendix B, 
table 9), the coefficient for SQBWAF is 0.129138 (nonsignificant), and the co- 
efficient for SQREGT is 0.21368 (significant). The effect of these two variables 
is to break the number of years since disturbance into two components: a 
significant increase when there is no budworm defoliation and a nonsignifi- 
cant increase when there is budworm. 

Figure 9 illustrates the concepts of optimum aspect and amplitude. Both 
are calculated from the interaction of slope and aspect using procedures devel- 
oped by Stage (1976). The mathematical determination of the highest and | 
lowest points along the curve indicates the general directions of the optimum 
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Figure 8—The effect of western spruce budworm and time on 
predicted probability of stocking for a Pseudotsuga menziesii/ 
Physocarpus malvaceus habitat type. 
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Figure 9—Illustration of optimum aspect and amplitude for the 
probability of stocking. 
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and poorest aspect for the dependent variable being predicted. The optimum 
aspect is the aspect with the highest point along the curve. In figure 9 it is 
17 degrees—indicating a general direction of north to northeast for the high- 
est probability of stocking for the Abies grandis/Clintonia uniflora habitat 
type. The poorest aspect is 180 degrees from the optimum aspect. 

Amplitude is the depth of the curve from the optimum aspect to the poor- 
est aspect. The larger the amplitude, the larger the difference between the 
optimum and poorest aspects. Amplitude can be calculated using any slope— 
we use 30 percent throughout this paper. The amplitude shown in figure 9 
is 0.34. This means the probability of stocking differs by 0.34 between the 
optimum and poorest aspects for a 30 percent slope. 

There is another important point about the calculation of amplitude when 
the form of the equation is 


P = (1+e28*))-1 


Amplitude is calculated with the curve centered on a probability of 0.5. 
Since probabilities are asymptotic to 0.0 and 1.0, amplitudes are at their 
maximum at a probability of 0.5. 


Analysis of data now excludes nonstocked plots. Given that a plot is stocked, 
we are interested in predicting how many trees are established, species compo- 
sition, and seedling heights. Distribution of the number of trees on stocked 
plots is shown in figure 10. The highest percentage is for plots that have 1 tree, 
followed in decreasing order by 2, 3, 4, and so on, up to 213 trees per plot. 
The result is a reversed J-shape distribution. 


Percent 


a 
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Number of Trees Per Stocked Plot 
Figure 10—Distribution of the number of regeneration-size trees 


per stocked plot. The tail of the distribution continues past 20 trees 
up to 213 trees on a plot. 
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The distribution of trees per plot is the important variable to reproduce, 
rather than the mean. Predictions of mean values with linear regression 
would not properly simulate the natural occurrence of regeneration because 
the categories with the highest probability of occurrence would seldom be 
predicted. Another method that simulates the distribution of the number of 
trees per plot is used. 

A two-step analysis was used to develop a Weibull cumulative density 
distribution equation. First, the data were categorized as follows: 


4 Habitat type series (Pseudotsuga menziesii, Abies grandis, 
Thuja/Tsuga, and Abies lasiocarpa) 


X 4 Aspects (north + 45 degrees, east, south, and west) 


X 3. Time periods (2-7 years since last disturbance, 8-12 years, 
13-20 years) 


X 2  Budworm defoliation histories (0-2 years, 23 years). 


Of the possible 96 categories, 73 had adequate data (225 plots). A Weibull 
equation was developed for each category, yielding 73 estimates of B and C for 
equation 2. Averages of variables by category were also calculated—average 
aspect, average slope, and so on. 

The second step in the analyses used linear regression to estimate Weibull 
parameters B and C as a function of site conditions. The resulting equations 
are shown in appendix B, table 10. Slope, aspect, elevation, habitat type series, 
number of years without budworm defoliation, and number of years with 
budworm defoliation are important predictors of the B parameter. For the 
C parameter, important predictors are elevation, number of years without 
budworm defoliation, and number of years with budworm defoliation. 


It is necessary to predict the number of species on stocked plots because all 
species do not occur on a given ¥0o-acre plot. This distribution has a reversed 
J-shape as shown in figure 11. The highest percentage is for 1 species, fol- 
lowed in decreasing order by 2, 3, 4, 5, 6, and 7 species. Seven species on a 
Yaoo-acre plot occurred only once. Separate equations were developed to pre- 
dict the probability of 1, 2, 3, 4, 5, and 6 species on stocked plots. Coefficients 
are given in appendix B, table 11. 

Prediction of the number of species on a stocked plot is conditional on the 
number of trees on the plot. For example, one tree on the plot could only 
result in one species on the plot; two trees, two or fewer species; three trees, 
three or fewer species; and so on. 


Advance regeneration is defined as best trees that germinated more than 
3 years prior to the year of harvest. Advance regeneration is commonly found 
beneath the canopy of mature stands in the Northern Rocky Mountains. Many 
of these trees survive harvest, site preparation, and the sudden exposure to 
more sunlight. Advance regeneration was included in the model as a sepa- 
rate class of regeneration so that future refinements to the Prognosis Model 
can account for growth response to release, important disease interactions 
with suppressed regeneration, and so on. Equations to predict the probabil- 
ity of advance regeneration were developed by species; coefficients are given 
in appendix B, table 12. 

The occurrence of advance regeneration is related to shade tolerance of 
the species. Shade-intolerant species have lower probabilities than shade- 
tolerant species. Occurrence is also related to the relationship of species to 
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Figure 11— Distribution of the number of species per stocked plot. 


habitat type. For example, Douglas-fir is considered shade tolerant on habi- 
tat types where it is the climax species, but is shade intolerant where it is 
not a climax species. 

Mechanical or burn site preparation generally decreases the probability 
of advance regeneration, with burn site preparation decreasing the probabil- 
ity more than mechanical. Because this study was retrospective, we could 
not determine whether advance regeneration was deliberately destroyed or 
whether none was present and site preparation was done to increase subse- 
quent regeneration success. 

The probability of advance regeneration decreases over time. This could 
be due to mortality over time or poor response to release. Trees that respond 
poorly to release could, with time, become subordinate to a subsequent tree 
of the same species. 


Trees up to 3 years old at the time of harvest were recorded as subsequent 
regeneration. This was for four reasons. First, there could be some error in 
aging small trees; accuracy was +1 year. Second, the year of harvest entered on 
historical records was usually the year of completion and, if it took more than 
1 year to harvest the stand, some areas would have been harvested before 
others. Third, a 3-year-old tree would not have suffered major physiological 
setbacks from suppression by overstory trees. Fourth, any trees planted im- 
mediately after harvest would appear to be advance regeneration. 

Equations for the probability of subsequent best trees are given in appen- 
dix B, table 18. The occurrence of subsequent regeneration differs from that 
for advance regeneration. Mechanical or burn site preparation increases 
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subsequent regeneration over untreated plots. In general, coefficients for 
mechanical and burn site preparation are also larger for shade-intolerant 
species versus shade-tolerant species. 

The number of years since last disturbance increases the probability of subse- 
quent species. The rate of increase for plots with spruce budworm defoliation is 
less than for plots without defoliation for three host species—Douglas-fir, 
grand fir, and subalpine fir. Spruce budworm defoliation in the 5 years prior 
to harvest also decreases the probability of subsequent grand fir. 

Residual overstory basal area decreases the probability of subsequent re- 
generation for every species except western hemlock. And if the species left 
in the overstory is the same one being predicted, the probability of subsequent 
regeneration is significantly increased for all species. 


Excess regeneration is all regeneration not chosen as best trees on the plot. 
The presence of an excess tree of a species is conditional on a best tree of the 
same species being established on the plot. Excess trees are distributed quite 
differently from best trees. Species that become established at high densities 
have a larger proportion of excess trees than do species that become estab- 
lished at low densities. For example, grand fir and western hemlock often 
overstock plots on the Tsuga heterophylla habitat type series. 

A consistently important independent variable increasing the probability of 
excess regeneration is the same species in the overstory (appendix B, table 14). 
The largest coefficient is for lodgepole pine and the smallest is for grand fir. 


Equations predicting heights of regeneration were developed for advance 
best trees, subsequent best trees, and excess trees. Various combinations 
and transformations of dependent and independent variables were tried. 
For advance and subsequent regeneration, the best dependent variable was 
always the natural log of tree height. Heights of excess trees are determined 
from Weibull distributions. 

Equations for heights of advance and subsequent best trees use tree age 
as an independent variable. This necessitates determination of tree age at 
the end of the Prognosis Model cycle. Weibull equations were developed that 
represent the distribution of the number of years from harvest to germina- 
tion of established seedlings. This is a negative value for advance regenera- 
tion and a positive value for subsequent regeneration. Once “delay to ger- 
mination” is known, tree age can be calculated. 


Heights of Advance Regeneration—Weibull distributions for choosing 
the delay to germination for advance best trees are given in appendix B, 
table 15. Delay to germination is dependent on species, residual overstory 
basal area, and spruce budworm defoliation during the 5 years before har- 
vest. Age is the number of years from germination to the end of the Prognosis 
Model cycle. Once tree age is determined, equations given in appendix B are 
used to predict tree heights. 

A total of 5,648 advance best trees were used to develop the height equa- 
tions in appendix B, table 16. Age is the single most important predictor of 
tree height. The coefficient for age is positive for all species; however, these 
coefficients are not as large as those for subsequent species. This is because 
the advance trees were suppressed prior to the harvest. 

Amplitude values for regeneration heights given in appendix B, tables 16 
and 18, are interpreted differently than the amplitude values in probability 
equations. For tree heights, amplitude is tree height at the optimum aspect 
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divided by tree height at the poorest aspect. For example, an amplitude of 
1.82 for advance lodgepole pine means that trees will be 1.82 times taller at 
the optimum aspect of 31 degrees than they are at the poorest aspect of 211 
degrees. 


Heights of Subsequent Regeneration—The calculation of heights for 
subsequent regeneration begins by determining the number of years from 
the last disturbance to seedling germination. Coefficients for these Weibull 
distributions are given in appendix B, table 17. Important variables are the 
number of years since last disturbance, spruce budworm defoliation after 
harvest, and species. 

Figure 12 shows cumulative Weibull distributions for the number of years 
of delay to seedling germination for plots 2-7, 8-12, and 13-20 years since 
disturbance. The species illustrated is grand fir with no budworm defoliation. 

Equations for predicting heights of subsequent best trees are given in appen- 
dix B, table 18. A total of 9,239 subsequent best trees were used in these 
analyses. As with advance regeneration, the most important independent 
variable is tree age. 

Spruce budworm defoliation history was not important for predicting height 
of host species. This finding agrees with Carlson (1988), who found that 
height growth of subsequent regeneration in even-age stands was not affected 
by spruce budworm and defoliation was light, even though adjacent overstory 
host trees were moderately to heavily defoliated. 

Increasing residual overstory density reduced heights for all species. Shade- 
intolerant species, such as ponderosa pine and lodgepole pine, had large 
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Figure 12—Weibull cumulative frequency distributions for delay to 
germination for subsequent grand fir regeneration. Curve A is for 
2-7 years time since disturbance, curve B for 8-12 years, and 
curve C for 13-20 years. 
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negative coefficients, while shade-tolerant species, such as western redcedar 
and western hemlock, had smaller negative coefficients. 

Site preparation increased heights for most of the shade-intolerant subse- 
quent species, but had no effect or a negative effect on shade-tolerant subse- 
quent species. 

Heights of subsequent regeneration increased as the number of trees on 
the plot increased. This seems counter-intuitive, because increasing compe- 
tition should decrease tree growth. However, increasing seedling density 
could indicate better microsites. 


Heights of Excess Regeneration—Heights of excess trees were not meas- 
ured in the field. However, a height must be assigned to all tree records being 
passed to the Prognosis Model. Two things are known about excess trees. 
First, an excess tree cannot be taller than the best tree of the same species 
on the plot. Second, an excess tree is at least as tall as the minimum estab- 
lishment height listed in table 1. 

Rather than assume a linear relationship between minimum establishment 
height and height of the tallest tree of the same species on the plot, Weibull 
distributions were developed for the heights of best trees. These distribu- 
tions were used to define the shape of the distribution for heights of excess 
trees. Coefficients are given in appendix B, table 19. Important variables 
for the distributions are species and years since last disturbance. 


HOW THE PROGNOSIS MODEL AND REGENERATION 
MODEL INTERACT 


The Prognosis Model is an individual-tree, distance-independent forest 
growth and yield model. Originally developed for northern Idaho and adjacent 
forests in Montana and Washington (Stage 1973; Wykoff and others 1982), the 
Prognosis Model has been calibrated for other forest regions. Geographic 
versions of the Prognosis Model are called variants. 

An inventory of trees from one or more plots represents a stand in the 
Prognosis Model. The inventory is maintained in a tree list where attributes 
of inventoried trees are kept—d.b.h., height, species, and so on. One of the 
attributes in the tree list is the number of trees per acre each tree represents. 
The trees per acre represented by a tree is initially a function of the inven- 
tory plot size and the number of plots sampled in the stand. Each tree sampled 
on a 1-acre plot represents 1.0 trees per acre, each tree sampled on two 1-acre 
plots represents 0.5 trees per acre, each tree sampled on ten ¥00-acre plots 
represents 30 trees per acre, and so on. 

The Prognosis Model projects the increment of individual trees over time. 
Time steps in the Prognosis Model are called cycles—cycles can vary in length, 
the default being 10 years. Prognosis Model harvests and thinnings, if any, 
are done at the beginning of the cycle; then growth and mortality are pre- 
dicted for each tree. Each tree receives a predicted increment in d.b.h. and 
height as well as a change in crown ratio. Mortality is simulated as a reduc- 
tion in the trees per acre that the tree represents in the stand. 

Plot attributes that change are updated at the end of the cycle for reporting 
purposes and for use by other submodels of the Prognosis Model; for example, 
plot overstory density and plot overstory species composition. The regen- 
eration model makes extensive use of plot information. 

Use of the regeneration model can be scheduled at the end of any Prognosis 
Model cycle. The regeneration model uses the updated inventory to predict 
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Step 1: Initiate 
Regeneration 
Model 


Step 2: Determine 
Plot Site 
Preparation 


regeneration that occurs during the cycle. The new trees are added to plots 
that correspond to the inventory plots in the Prognosis Model. The use of 
plot information provides some spatial resolution to the regeneration model 
because each plot can have unique attributes such as slope, aspect, habitat 
type, site preparation, overstory density, and overstory species composition. 

New regeneration is added to the Prognosis Model tree list at the end of the 
cycle in what is called a regeneration tally. A tally of regeneration is similar 
to the results of taking a regeneration survey for a stand; the regeneration 
process continues for a number of years before it is quantified by a regeneration 
survey. 

The regeneration period is the number of years during which regeneration 
becomes established as a result of a disturbance. For modeling purposes, the 
regeneration period is set at 20 years. Any time there is a Prognosis Model 
cycle boundary during the 20 years of a regeneration period, the regenera- 
tion model will be invoked to predict new regeneration. Multiple tallies during 
these 20 years are called a tally sequence. 

Regeneration that occurs after the regeneration period is called ingrowth. 
Ingrowth is both the result of succession by shade-tolerant species and the 
continued regeneration of trees into gaps in the tree canopy. 


STEPS IN THE REGENERATION MODEL 


The computer steps for the regeneration model are depicted in figure 13. 
The computer code is written in FORTRAN. Each step is briefly described 
in this section. 


The regeneration model is invoked at the end of a Prognosis Model cycle. 
At this point in the Prognosis Model, silvicultural prescriptions for the cycle 
have been implemented, growth has been predicted, mortality has been ac- 
counted for, and plot/stand statistics have been updated. 

Before the regeneration model is invoked, the tree list is checked to ensure 
the tree list has enough space to hold new records. Space is provided by 
compressing the tree list. This feature of the Prognosis Model combines 
similar tree records as discussed by Wykoff (1986). 

In the regeneration model, computer storage is allocated and prepared for 
use. Plot replication, if needed, is done at this time. Plot replication ensures 
there are an adequate number of plots (default = 50) to dampen the stochas- 
tic effects of any one plot on the projection. The effect is to predict regenera- 
tion closer to the mean expectation. Fifty plots are enough to average out 
the effects of unusual events, yet is still efficient use of computer time. Model 
users can change the default minimum number of plots to as few as 20. 


Each plot is assigned a site preparation—mechanical, burn, or untreated. 
The percentages of site preparations are either specified by the user or are 
determined from default equations derived from the data. User-specified 
percentages take precedence over default equations. Once percentages are 
determined, site preparations are assigned at random to the plots. 

Users have two ways to assign site preparations to specific plots. Site prepa- 
ration can be coded on the inventory records, or it can be specified on special 
computer records at the time the Prognosis Model is initiated (see Ferguson 
and Crookston 1991). A category for roads is included. Road site prepara- 
tion is all road cuts, all road fills, and unmaintained roadbeds. | 
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No 
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Figure 13—Steps in the Regeneration 
Establishment Model. 
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Step 3: Compute 
Years Since Last 
Disturbance 


Step 4: Calculate 
the Increment in 
Stocking 


Step 5: Estimate. 
Number of Trees 


For each plot, the number of years since last disturbance is computed—the 
result is dependent on the method of site preparation. Untreated plots were 
last disturbed at the time of harvest. Mechanically disturbed plots were last 
disturbed at the time of mechanical site preparation, and burned plots were 
last disturbed at the time of burning. 


The increment in stocking is the increase in the probability of stocking from 
one time period to another. Figure 14 illustrates the increment in the prob- 
ability of stocking for successive time periods of 5, 10, 15, and 20 years. The 
number and length of time periods can be altered by controlling cycle lengths 
in the Prognosis Model. Whenever there is a cycle boundary in the Progno- 
sis Model and a disturbance has occurred in the previous 20 years, a regen- 
eration tally will be predicted by the regeneration model. 

A calibration feature in the regeneration model is invoked when the stand 
inventory is taken within 20 years of a requested tally of regeneration. Calibra- 
tion adjusts the probability of stocking curve to coincide with the actual inven- 
tory at the time the inventory was taken. Then the increment in stocking from 
the time of inventory to the end of the cycle boundary is calculated. Details 
of the calibration feature are explained in a later section. 


The number of trees per plot is conditional on the plot being stocked. In the 
first tally of a tally sequence, the number of trees on the plot is determined 
from Weibull cumulative distribution functions developed from the data. A 
uniformly distributed pseudo-random number is used to choose the actual 
number of trees from Weibull distributions. The number of trees is rounded 
to the nearest integer. 
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Figure 14—Increments in probability of stocking at 5, 10, 15, and 
20 years time since disturbance. | 
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Step 6: Predict 
Number of 
Species 


Step 7: Calculate 
Probability of 
Species 
Occurrence 


Step 8: Estimate 
Heights 


Step 9: Add 
Planted Trees 


For other tallies of a tally sequence, only the incremental number of trees 
becoming established on the plot during that time period is added to the tree 
list. 


The number of species is also conditional on the plot being stocked. Logistic 
regression equations predict the probability of 1, 2, 3, 4, 5, and 6 species on 
the plot. The number of species cannot exceed the number of trees chosen 
in step 5. 

A uniformly distributed pseudo-random number is used to make a discrete, 
but unbiased, choice of the number of species on the plot. The probabilities 
for each number of species are calculated and totaled. The total is divided 
back into the probabilities so that the sum of the probabilities equals 1.0. The 
adjusted probabilities are accumulated within the interval [0,1], and a uni- 
formly distributed pseudo-random number is compared to the accumulated 
probabilities. The number of species on the plot is the one for which the accu- 
mulated probability first exceeds the pseudo-random number. 


There are 10 species in the regeneration model and each can occur in any 
of three ways: as advance best trees, as subsequent best trees, and as excess 
regeneration. Advance trees germinated more than 8 years prior to the har- 
vest date. Subsequent trees germinated after the cutoff date for advance 
trees. Excess trees were established on the plot, but were not aged. 

First, the probabilities of advance and subsequent best trees are predicted 
by species (probabilities can be 0.0 when certain species do not occur on some 
habitat types or in some geographic locations). The advance and subsequent 
distributions are accumulated within the interval [0,1], and pseudo-random 
numbers are used to choose the same number of species on the plot as was 
chosen in step 6. 

Next, the probability of excess regeneration is predicted by species. Species 
of remaining trees on the plot are determined from the excess distribution. 
Probabilities in the excess regeneration distribution can differ considerably 
from the combined advance and subsequent distribution. 


Heights of best trees are predicted from linear regression equations. Separate 
equations were developed by species for advance and subsequent best trees. 

Heights are predicted for each best tree on the plot. Tree heights are var- 
ied by adding or subtracting a random proportion of the standard error of 
the estimate to predicted tree height. A pseudo-random number determines 
whether to add or subtract. A second pseudo-random number in the inter- 
val [0,1] determines the proportion of the standard error to add or subtract. 

Excess tree heights are determined from Weibull cumulative distribution 
functions that were developed from the heights of best trees. A pseudo- 
random number is used to assign heights for excess trees between the mini- 
mum establishment height and the height of the best tree of the same species 
on the plot. 


The regeneration model predicts natural regeneration. Planted trees are 
an addition to predicted natural regeneration. Model users specify species, 
density, year of planting, and survival of planted trees. The default option 
is for trees to be planted uniformly throughout the stand, but options allow 
for planted trees to vary in density in relation to residual overstory basal 
area (see Ferguson and Crookston 1991). 
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Step 10: Identify 
Best Trees 


Step 11: 
Accumulate Stand 
Statistics 


Step 12: Pass Tree 
Records 


Step 13: Print 
Regeneration 
Summary 


Heights of planted trees are calculated using equations developed for 
heights of subsequent trees. Age of planted trees is the number of years 
from planting to the end of the Prognosis Model cycle plus the age of the 
tree when planted. 


Best trees are chosen from those that are established on the plot using the 
rules described in the section on field procedures. Best trees are assigned a 
lower priority for removal in thinning operations simulated by the Prognosis 
Model (see Wykoff and others 1982, pp. 16-17). Excess trees have a higher 
priority for removal. 

Planted trees are included in the selection of best trees. Planted trees could 
be taller than natural regeneration of the same species; however, delays in 
planting could result in the natural regeneration being taller than planted 
trees, so planted trees are not always chosen as best trees. 


This is primarily a bookkeeping step for the purpose of later printing the 
regeneration summary. Predictions by plot are stored for averaging that 
produces stand statistics. 


Tree records are created during this step. Tree.records created by the re- 
generation model are added to the Prognosis Model tree list. This step uses 
both the probability of stocking and the results of stocked plot predictions to 
assign a “trees-per-acre” value to new tree records. The trees-per-acre value 
assigned to each tree record is 


TPA = (PS * TPSP * 300)/N (4) 
where 


TPA = trees per acre 
PS = plot probability of stocking 
TPSP = trees per stocked plot 
300 = inverse of plot size used in sampling 
N = number of plots and replicates being processed. 


Each best tree on the plot becomes a separate tree record, so TPSP is 1.0 
for best trees. Excess trees are accumulated by species, then tree records 
are created so that TPSP is <5 trees. 

The value of PS for planted trees is not the probability of stocking; rather, 
it is the density of planting specified by the user. Usually one tree record 
for each species being planted is created for each plot. Multiple records are 
created when high densities are planted (a tree record represents >10 trees 
per acre). 


A regeneration summary is printed each time the regeneration model is 
invoked. The first tally in a tally sequence includes a heading that shows 
the stand identification, site preparation by year performed, and number of 
plots by habitat type group. 

Each tally shows the cumulative probability of stocking and three summaries 
of regeneration by species. The first summary is for all trees regenerating 
during the cycle. This summary is similar to a conventional stand exam. The 
second summary lists best trees regenerating during the cycle. Densities 
and average heights are listed by species for best trees. Focusing on best 
trees gives a more realistic picture of potential crop trees in the stand. For 
example, one could determine the number of Douglas-fir that could be favored 
in the stand through a thinning operation. 
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Example 
Projection 


Comparison of 
Habitat Types 


The third summary shows all trees less than 3.0 inches d.b.h. that are in 
the Prognosis Model tree list. This provides a view of all regeneration cur- 
rently in the stand. Trees may be from a stand inventory, from a previous 
tally of the regeneration model, and from the current tally of the regenera- 
tion model. 


MODEL BEHAVIOR 


In this section we use the regeneration model to show example projections. 
Examples include the default conditions that have been used to illustrate 
equations, a comparison of five habitat types, and an example of the effects 
of western spruce budworm. 


This projection is the same stand used to illustrate equations—a hypothetical 
clearcut in the Idaho Panhandle National Forests, Abies grandis/Clintonia 
uniflora habitat type, north aspect, 30 percent slope, 3,500 feet elevation, no 
planting, and no site preparation. 

The computer file needed to make this projection (called keyword records 
in the Prognosis Model system) is shown in figure 15. Refer to Wykoff and 
others (1982) and Ferguson and Crookston (1991) for coding instructions. The 
default cycle length for the Prognosis Model is 10 years, so regeneration will 
be predicted at 10 and 20 growing seasons since the disturbance in 1990. 

Output from the regeneration model is shown in figure 16. New tree records 
are added at the end of the Prognosis Model cycle. There have been 10 growing 
seasons between the “spring” of 1990 and the “fall” of 1999. The new regenera- 
tion is added at the end of the cycle so these trees are part of the stand when 
silvicultural prescriptions are simulated at the beginning of the next cycle. 

After 10 growing seasons, the probability of stocking is 0.7380 and regenera- 
tion density is 772 trees per acre. Of the 772 trees, 562 are best trees. Spe- 
cies composition of best trees is 57 percent grand fir, 22 percent Douglas-fir, 
9 percent western white pine, 5 percent western larch, and less than 5 percent 
each of Engelmann spruce, lodgepole pine, and ponderosa pine. 

After 20 growing seasons, the probability of stocking is 0.8631 and 309 more 
trees have become established, 209 of which are best trees. There are now 
989 trees <3 inches d.b.h. in the Prognosis Model tree list. Some of these 
trees are from tally 2 and some are from tally 1. 


Here we show how regeneration model projections vary with habitat type 
and elevation. Five projections were made, changing only habitat type and 
elevation between each projection. Other variables held constant were Idaho 
Panhandle National Forests, north aspect, 30 percent slope, no planting, no 
residual overstory basal area, no site preparation, and 10 years since harvest. 
The five projections were a Pseudotsuga menziesii/Physocarpos malvaceus 
habitat type at 3,000 feet elevation, an Abies grandis/Clintonia uniflora habitat 
type at 3,500 feet, a Tsuga heterophylla/Clintonia uniflora habitat type at 
4,000 feet, an Abies lasiocarpa/Clintonia uniflora habitat type at 5,000 feet, 
and an Abies lasiocarpa/Xerophyllum tenax habitat type at 5,500 feet. 

Results of these projections are shown in table 6. The Pseudotsuga menziesii/ 
Physocarpos malvaceus habitat type has the lowest probability of stocking 
(0.2596) and density (195 trees per acre) of the five projections. The number 
of species present in the stand is limited because of the habitat type. The 
Abies grandis/Clintonia uniflora habitat type has a probability of stocking 
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SCREEN 


STDIDENT 
ABGRCLUN EXAMPLE KEYWORDS FOR REGENERATION ESTABLISHMENT MODEL 
DESIGN 10.0 0.0 
STDINFO 104.0 520.0 0.0 0.0 30.0 35/0 
INVYEAR 1990.0 
NUMCYCLE 350) 
NOTREES 
ESTAB 1990.0 
END 
PROCESS 
STOP 
Figure 15—Example keyword records for projecting a stand in the 
Prognosis Model. Results of the regeneration model for this projec- 
tion are shown in figure 16. 
STAND GROWTH PROGNOSIS SYSTEM VERSION 6.10 -- INLAND EMPIRE 
REGENERATION ESTABLISHMENT MODEL VERSION 2.0 
STAND ID: ABGRCLUN MANAGEMENT CODE: NONE 
SITE PREP SUMMARY PLOT HABITAT TYPE SUMMARY 
PREP:NONE MECH BURN SERIES:<--- DF --> <--- GF --> C WH <--- AF & MH ---> 
YEAR:1990 1990 1990 GROUP?) 12083) “4° 5) 67) (ss 9) 10) Aa 2 Sh aaa Seale 


PCT: 68 24 8 Peon (0) Ol OF Oh Oi Oslo OO @ WO © © @ @"@ 


TALLY 1 AT 10 YEARS. PROBABILITY OF STOCKING IS 0.7380 IN THE FALL OF 1999. 


SUMMARY OF ALL 


TREES REGENERATING 
DURING THIS TALLY. 


SUMMARY OF BEST 


TREES REGENERATING 
DURING THIS TALLY. 


TREES 


< 3.0 IN. DBH 


BEING PROJECTED BY 
THE PROGNOSIS MODEL 


TREES % OF TREES % OF AVERAGE TREES % OF 
SPECIES /ACRE TOTAL /ACRE TOTAL HEIGHT /ACRE TOTAL SPECIES 
WP 49. 6. 49. 9. 4.5 49. 6. WP 
L 30. 4. 30. 5. Lege 30. 4. L 
DF 193. 25. 124. 22. 4.7 193. 25. DF 
GF 461. 60. 319%. 57. 3.6 461. 60. GF 
WH O. 0. O. 0. 0.0 QO. OF WH 
Cc QO. QO. QO. QO. 0.0 QO. QO. Cc 
LP 9. 1. Oe 2. 2.3 9. 1. LP 
Ss 22. 3. 22. 4. 2.6 22. 3. s 
AF O. QO. QO. QO. 0.0 QO. QO. AF 
PP 9. 1. 9. 2. 3.5 9. 1. PP 
a QO. QO. QO. QO. 0.0 QO. O. SS 

772. 562. 772. 


TALLY 2 AT 20 YEARS. 


SUMMARY OF ALL SUMMARY OF BEST TREES < 3.0 IN. DBH 
TREES REGENERATING TREES REGENERATING BEING PROJECTED BY 
DURING THIS TALLY. DURING THIS TALLY. THE PROGNOSIS MODEL 
TREES OF TREES @ OF AVERAGE TREES % OF 
SPECIES /ACRE TOTAL /ACRE TOTAL HEIGHT /ACRE TOTAL SPECIES 
WP 24. 8. 24. 11. 3.3 69. Uo WP 
L 6. 2. 6. 3. 4.2 33. 3. L 
DF 56. 18. 40. 19. 3.3 220. 22. DF 
GF 206. 67. 124. 59. 2.2 618. 62. GF 
WH QO. QO. 0. QO. 0.0 QO. QO. WH 
Cc QO. QO. O. QO. 0.0 QO. QO. Cc 
LP QO. QO. 0. QO. 0.0 9. 1. LP 
Ss 14. 5. 14. 7. 2.0 36. 4. BS) 
AF 1. QO. abe QO. Oronl, 1. QO. AF 
PP 1. QO. 0G QO. 4.7 3. QO. PP 
=m QO. QO. 0. QO. 0.0 O. QO. os 
309. 209. 989. 


PROBABILITY OF STOCKING IS 0.8631 IN 


Figure 16—Regeneration model summary table. Keywords shown in figure 
15 were used to make this projection. 
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THE FALL OF 2009. 


Example of 
Spruce Budworm 
Effects 


Table 6—Predicted probability of stocking and trees per acre by species for five projections 
where only habitat type and elevation were changed. Variables held constant were 
the Idaho Panhandle National Forests, north aspect, 30 percent slope, no planting, 
no residual overstory basal area, no site preparation, and 10 years since harvest 


Habitat type and elevation 
PSME/PHMA' ABGR/CLUN TSHE/CLUN ABLA/CLUN' ABLA/XETE 


Stocking 3,000 ft 3,500 ft 4,000 ft 5,000 ft 5,500 ft 
Probability 
of stocking 0.2596 0.7380 0.9015 0.6272 0.6671 
Trees/acre 195 772 2,248 746 768 
Trees/acre 
by species 
WP? 0 49 129 16 22 
le 18 30 5 18 82 
DF 148 193 222 96 225 
GF 0 460 801 236 61 
WH 0 0 682 0 0 
C 0 0 268 0 0) 
LP 10 9 16 3 Al 
iS) 0 22 5 50 79 
AF 0 0 120 327 228 
PP 19 9 0 0 0 


'See table 2, footnote 1 for habitat type abbreviations. 
See table 1 for species abbreviations. 


of 0.7380 and 772 trees per acre. Species composition is more diverse with 
grand fir having the highest density. 

The Tsuga heterophylla/Clintonia uniflora habitat type has the highest prob- 
ability of stocking (0.9015) and density (2,248 trees per acre) of the five habi- 
tat types projected. The Abies lasiocarpa/Clintonia uniflora habitat type regen- 
erates well to subalpine fir (327 trees per acre). The probability of stocking on 
this habitat type is 0.6272, and density is 746 trees per acre. 

On the Abies lasiocarpa/Xerophyllum tenax habitat type, the probability of 
stocking and density are close to the Abies lasiocarpa/Clintonia uniflora habitat 
type. The probability of stocking is 0.6671, and density is 768 trees per acre. 
The Abies lasiocarpa/Xerophyllum tenax habitat type has more western larch, 
Douglas-fir, and lodgepole pine. Grand fir is less frequent. 

These five examples illustrate the diversity of results that can be expected 
among habitat types. Growth and mortality predictions in the Prognosis 
Model are also sensitive to habitat type, which helps represent the diversity 
of ecosystems found in the Northern Rocky Mountains of the Western United 
States. 


The effects of western spruce budworm on regeneration success have not 
been dramatic in the separate equations shown thus far. However, the cumula- 
tive effects of budworm can be quite important, as will be illustrated by exer- 
cising the model with and without budworm. 

Two hypothetical stands were chosen for illustration. The first is a clearcut 
in the Lolo National Forest in Montana at 4,500 feet elevation, north aspect, 
30 percent slope, no planting, no site preparation, Pseudotsuga menziesii/ 
Symphoricarpos albus habitat type, and 10 years since harvest. The second 
stand is the same except the elevation is 5,500 feet and the habitat type is 
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Abies lasiocarpa/Carex geyeri. These two stands were projected with and 
without budworm defoliation for the time period 5 years before harvest 
through 10 years after harvest. 

Results of these four projections are shown in table 7. Budworm decreases 
the probability of stocking from 0.4109 to 0.3485 for the Pseudotsuga menziesii/ 
Symphoricarpos albus habitat type and from 0.4391 to 0.3005 for the Abies 
lasiocarpa/Carex geyeri habitat type. Density also decreases in the defoliated 
stands. The Pseudotsuga menziesii/Symphoricarpos albus habitat type has 
310 trees per acre without budworm compared to 250 with budworm. The 
Abies lasiocarpa/Carex geyeri habitat type has 450 trees per acre without 
budworm compared to 295 with budworm. The effect of budworm defoliation 
on stocking and density can be to reduce stocking below acceptable levels. 

In addition, budworm defoliation changes species composition. The density 
of host species is reduced; for example, the density of Douglas-fir is 161 trees 
per acre without budworm compared to 123 with budworm in the Pseudotsuga 
menziesit/Symphoricarpos albus habitat type. In addition to changes in species 
composition, there can be a shift in the proportion of advance and subsequent 
regeneration within a species. 

Seedling heights are also affected by budworm. Host species are shorter. 
This is caused by a longer delay to germination of an established seedling or 


Table 7—Results of projecting two stands with and without budworm defoliation for the 
time period 5 years before harvest through 10 years after harvest. Results are 
shown as of 10 years 


Habitat type and elevation 


PSME/SYAL,' 4,500 ft ABLA/CAGE, 5,500 ft 
Stocking Budworm No budworm Budworm No budworm 
Probability 
of stocking 0.3485 0.4109 0.3005 0.4391 
Trees/acre 250 310 295 450 
Trees/acre 
by species 
WP? 0 0 0 0 
L 62 15 0 0 
DF 123 161 39 63 
GF 0 0 11 25 
WH 0 0 0 0 
Cc 0 0 0 0 
LP 29 32 105 164 
S 0 0 50 37 
AF 0 0) 90 161 
PP 36 42 0) 0 
Average heights 
of best trees (ft) 
L: 4.2 4.4 — — 
DF 4.9 4.2 4.0 3.9 
GF — — 2.8 3.0 
LP Sh7/ 3.4 Sal 3.0 
S) — — 25 2, 
AF “= — 4.1 3.8 
PP 5.7 5.6 _ = 


'See table 2, footnote 1, for habitat type abbreviations. 
See table 1 for species abbreviations. 
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Calibration 


by slower growth. Yet, notice that the average height of best trees for Douglas- 
fir in the Pseudotsuga menziesii/Symphoricarpos albus habitat type is 4.9 feet 
for the budworm projection and 4.2 feet without budworm (table 7). This 
occurs because budworm changes the proportion of advance and subsequent 
regeneration within a species. Budworm reduces the probability of subse- 
quent Douglas-fir (appendix B, table 13) but does not change the probability 
of advance Douglas-fir (appendix B, table 12). Thus, there is a greater pro- 
portion of taller advance Douglas-fir, increasing the average height for that 
species. 

Nonhost species can be taller as well. This effect could be due to earlier estab- 
lishment of nonhost species or better growth rates because of less competition. 

In summary, there are four effects of western spruce budworm represented 
in the regeneration model. First is a reduction in the probability of stocking. 
Second is reduced number of trees per stocked plot. The significance of reduc- 
tions in the probability of stocking (PS) and the number of trees per stocked 
plot (TPSP) becomes clear when densities are assigned using equation 4: 


TPA = (PS * TPSP * 300)//N 


PS and TPSP interact multiplicatively, so small reductions in either (or both) 
can have large impacts. 

The third effect of budworm is to change species composition of the new 
stand. Fewer host species regenerate in the new stand. The fourth effect 
is to change heights of regeneration. Nonhost species may become taller 
relative to host species that are under attack by budworm. Subsequent host 
species in defoliated stands may be shorter, primarily due to delays in be- 
coming established. However, average heights of host species may be taller 
with budworm defoliation than without defoliation. This happens when the 
proportion of advance regeneration for a species increases relative to subse- 
quent regeneration. The taller advance regeneration increases the average 
height for the species. 


OTHER MODEL FEATURES 


Most of the details about the regeneration model have been discussed. 
A few additional topics may be of interest to others developing a similar model. 


Model calibration is a process of adjusting the probability of stocking curve 
to coincide with an inventory. This adjustment takes place whenever the 
year of inventory is within the years the regeneration model is making pre- 
dictions. A stand more densely stocked than expected has less room for ad- 
ditional regeneration, while a stand less densely stocked than expected has 
more room for regeneration. 

We considered using the ratio of stocked plots to calibrate the probability 
of stocking. However, this calibration method would not work for inventories 
with plot sizes other than Y00-acre. 

This drawback was eliminated by developing an equation predicting the 
proportion of plots that are stocked from the number of trees per acre. The 
equation is 


RATIO = [1 + exp — (-5.17397+0.85131*In(TPA)) |" (5) 
where 


RATIO = proportion of plots that are stocked 
TPA = trees per acre in the stand <3 inches d.b.h. 


33 


Consistency of 
Predictions 


This equation has an r-square of 0.7694 and was developed from 362 stands 
in which 10 or more ¥o0-acre plots were sampled and on which all established 
regeneration was recorded. The relationship between stocking and density 
is shown in figure 17. 

The predicted ratio of stocked plots is then used to mathematically adjust 
the probability of stocking equations given in appendix B, table 9. The ad- 
justment is made to 2B,X; in the probability equation 


P = (1+e-8%))1 


so that the probability remains bounded in the interval [0,1]. 

Figure 18 illustrates the calibration procedure. Line A is the probability 
of stocking curve generated from the Abies grandis series equation in ap- 
pendix B, table 9. Line B represents the same stand that was inventoried 
at 3 years since harvest and is more densely stocked than expected. Line C 
represents the same stand as line A, but it is less densely stocked than ex- 
pected at 3 years. Calibration adjusts lines B and C to coincide with the in- 
ventories (points b and c). Adjusted predictions (lines B and C) are still 
bounded within the interval [0,1]. 

A second part of calibration concerns creating new tree records. If the ac- 
tual number of trees on the plot exceeds the predicted number, no new tree 
records are created for that plot. 


The regeneration model uses pseudo-random numbers to make unbiased, 
discrete choices from a number of possibilities; for example, the species that 
will occupy a plot. This feature mimics stochastic variation seen in nature. 
Regeneration in a stand would vary depending on the year it was harvested 
because of differences in weather, seed crops, and animal/insect activity 
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Figure 17—The relationship of trees per acre to stocked plot ratio. 
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Figure 18—Use of the calibration feature in the Regeneration Es- 
tablishment Model. Line A is the predicted probability of stocking 
curve. Point b represents an inventory from the same stand but 
stocking is better than predicted. This information is used to ad- 
just the probability stocking curve upward to line B. Point c repre- 
sents an inventory from the same stand as line A but stocking is 
poorer than predicted. The probability of stocking curve is ad- 
justed downward (line C). 


from year to year. Variation can also be expected among inventories of the 
same stand. If different starting points are chosen for a transect or sampling 
grid, one would expect different results from the inventory. 

The use of pseudo-random numbers in the regeneration model also resulted 
in inconsistencies among projections of the same stand where projections 
differed only in cycle lengths. For example, when two 5-year cycles are used 
instead of one 10-year cycle, the pseudo-random numbers are used in a differ- 
ent sequence between projections. Consistency was obtained by using pseudo- 
random numbers in the same order for all tallies within a tally sequence. 

To further improve consistency, species composition is held constant at 
10 years since last disturbance. This is necessary because of the relative 
increases and decreases over time in the probability of species presented in 
appendix B. Species composition is relatively stable at 10 years, and none 
of the species in the model has long delays to establishment. 

An example of consistency among runs is given in table 8. Three projec- 
tions are compared—four 5-year projections, two 10-year projections, and one 
20-year projection. The probability of stocking, density, and species composi- 
tion are nearly identical, regardless of cycle length. The small differences 
that do occur are attributable to very slight differences in probabilities when 
regeneration from one tally adds basal area to the stand in the next tally— 
the probability of stocking and species probabilities are shifted slightly. 


35 


Ingrowth 


Automatic Tallies 


Table 8—Example projections to compare consistency of various cycle lengths using the same 
stand. Projection A is four 5-year cycles, projection B is two 10-year cycles, and pro- 
jection C is one 20-year cycle 


Projection A Projection B Projection C 
Stocking Syr 10yr 15yr 20 yr 10 yr 20 yr 20 yr 
Cumulative 
probability 
of stocking 0.6186 0.7380 0.8126 0.8669 0.7380 0.8631 0.8601 
Cumulative 
trees/acre 582 772 943 1,086 772 1,081 1,076 
Cumulative 
trees/acre 
by species 
WP' 41 49 59 74 49 73 72 
(le 21 30 34 37 30 36 41 
DF 157 192 222 248 193 249 262 
GF 330 462 575 667 460 667 635 
WH 0 0 0 0 0 0 0 
C 0 0 0 0 0 0 0 
LP 7 8 14 20 9 9 20 
Ss 18 22 29 30 22 36 36 
AF 0 0 0 0 0 1 0 
PP 8 9 10 10 9 10 10 


'See table 1 for species abbreviations. 


Ingrowth is regeneration that occurs after the 20-year regeneration period. 
It is the continued succession by shade-tolerant species, and the regeneration 
of trees into gaps in the overstory canopy. An ingrowth tally of regeneration 
is predicted every 20 years if the regeneration model has not been invoked 
during that time, and no regeneration activities are scheduled for the next 
Prognosis Model cycle. 

Data collected for the regeneration model did not specifically address in- 
growth into undisturbed stands. However, the data contain many plots repre- 
senting conditions where ingrowth occurs. These are areas with fairly dense 
overstory canopies where no harvesting or site preparation had taken place 
(see table 5). The systematic arrangement of transects or sampling grids in 
selection, sanitation, or salvage harvests resulted in many of these plots being 
included in the sample. 

When predicting ingrowth, all plots are simulated as receiving no site prepa- 
ration. Ingrowth tree records are created and passed to the Prognosis Model. 
An assortment of tree sizes is added to the stand because trees up to 20 years 
old can be established. The trees then survive and grow based on tree, site, 
and stand characteristics. 

Simulation of ingrowth creates more realistic projections of stand develop- 
ment over time. It also allows projections that simulate succession over long 
time periods. The Prognosis Model with the regeneration model can be used 
in areas of forest management that require evaluation of uneven-age and 
old-growth silviculture. 


The Prognosis Model includes a feature that invokes an automatic scheduling 


of the regeneration model when trees are removed. Thinnings scheduled by the 
user may have opened the stand enough for regeneration to become established. 
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These situations are detected and automatic tallies scheduled when certain 
threshold values are exceeded. 

If 10 to 30 percent of either the trees per acre or the total cubic-foot vol- 
ume is removed, a single regeneration tally is scheduled. A tally sequence 
is scheduled when more than 30 percent of either the trees per acre or the 
cubic-foot volume is removed. Nothing is scheduled if removals are less than 
10 percent. The default values of 10 and 30 percent can be changed by the 
user. 


Site Preparation Site preparation is a very important predictor of regeneration success, 
species composition, and growth rates. Most model users either specify the 
percentages of site preparation in the prescription, or they code plot-by-plot 
site preparations on stand examination records. For model users who do 
not specify site preparation(s), equations were developed to predict default 
probabilities. Plots used to develop these equations came from stands not 
deliberately prepared mechanically or by burning. 

Default site preparation equations are given in appendix B, table 20. Only 
aspect, slope, residual basal area, elevation, and habitat type series are im- 
portant independent variables. 


Regeneration The Prognosis Model was first developed for northern Idaho and adjacent 

Model Variants portions of Montana and Washington (Stage 1973; Wykoff and others 1982). 
The Prognosis Model has also been calibrated for other geographic areas. 
Geographic versions are called variants. 

Prognosis Model variants that do not have a calibrated regeneration model 
have an abbreviated version of the regeneration model. In these variants, 
creation of tree records representing natural regeneration is specified by the 
user in a manner similar to planting (Ferguson and Crookston 1991). Users 
describe expected regeneration in terms of species and density. Expected 
regeneration is estimated by the user using historical data or experience. 


Regeneration Some Prognosis Model variants represent species that sprout from roots 

Sprouts or stumps of harvested trees. For these variants, subroutines that create 
tree records for sprouts were added to the regeneration model. Information 
about trees removed at the beginning of the cycle is stored for use at the end of 
the cycle when the number and size of sprouts is predicted. Each variant 
will be different because of the species and environment being represented, 
but once equations predicting the number and size of sprouts are developed, 
they can readily be inserted in the regeneration model. 


DISCUSSION 


Version 2 of the Regeneration Establishment Model became available 
with the release of version 6.1 of the Prognosis Model. Version 2 of the re- 
generation model covers five habitat type series in Montana and in central 
and northern Idaho. Abbreviated regeneration models are available with 
Prognosis Model variants covering other locations in the Western United 
States. 

The regeneration model is based on a sample size of 12,128 Yo0-acre plots 
from 537 stands. Data represent all aspects, 0-110 percent slopes, 0-390 
square feet of residual overstory basal area, 2,400-7,400 feet elevation, four 
site preparation methods, and 0-16 years of budworm defoliation. Stands 
were chosen using stratified random sampling procedures. Stratification 
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assured replicated coverage of important variables—habitat type series, site 
preparation, overstory density, and level of spruce budworm defoliation. 
Stands were randomly selected to ensure an unbiased choice among all can- 
didate stands. Plots within stands were made as independent as possible 
by recording plot-specific variables whenever feasible and by spacing plots 
throughout the stand. 

We feel that unbiased stand selection procedures are crucial to developing 
a regeneration model because results are very dependent on probabilities. 
For example, the probability of stocking depends on the proportions of stocked 
and nonstocked plots. These proportions could be easily changed if study 
areas are not chosen in an unbiased manner. 

The importance of random stand selection does not mean that case history 
studies of individual or small groups of stands are not needed. These studies 
provide very useful information on silvics of species, regeneration success, 
processes leading to successful regeneration, species composition, and so on. 
Case history studies helped us design this study, but it is inappropriate to 
use case history studies to develop a regeneration model because of potential 
problems with low number of stands in the sample, narrow range of years 
of treatment, and the possibility of unintended stand selection bias. 

Stands used to develop the regeneration model were chosen from a list of 
stands harvested by commonly used methods in the Northern Rocky Moun- 
tains. These stands were harvested as part of commercial timber removals 
on Ranger Districts and State and corporate ownerships. Thus, the regen- 
eration model should reflect what can be expected when actual silvicultural 
prescriptions are implemented. 

Retrospective examination procedures may not be suitable for developing 
regeneration models in all cases. Sometimes stand conditions prior to harvest 
may be important in determining regeneration in the next rotation. For exam- 
ple, the presence of serotinous lodgepole pine cones is important for lodgepole 
pine regeneration. 

Retrospective examination procedures used in this study have the drawback 
of sampling stands that are up to 20 years old. During that time, technology 
could improve regeneration success. The same concern is true of a prospec- 
tive study. Retrospective examinations can be conducted in a much shorter 
time frame—say 2 to 3 years to collect a sufficient amount of data. Retrospec- 
tive studies, however, can only be done if adequate historical records have 
been kept for each stand. 

Another important part of data collection is to sample from a range of treat- 
ment years. This spreads out the sporadic effects of variables not explicitly 
accounted for in the model—weather, seed crops, animals, insects, disease, 
and so on. Data collected over a range of treatment years should result in 
model predictions that are closer to long-term averages. 

An ecologically based land classification system was very important for 
model development. The regeneration model makes extensive use of habitat 
types and habitat type groups. Other alternatives would not have worked 
as well. For example, site index does not provide the necessary information 
to model regeneration success. Site index alone is not useful for determin- 
ing the species that can occupy the site. 

Ecological land classifications provide a logical framework within which 
to model regeneration success. The successional pathways to a climax for- 
est can be modeled using an ecological classification system. Management 
implications can be cataloged using the same system. The analyses and model 
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presented in this paper are ecologically based and should give reliable esti- 
mates of future stocking under various management strategies. 

Although we classified stands by regeneration method—clearcut, seedtree, 
shelterwood, and selection—many stands did not fit the classical definition, 
and there was usually substantial variation within stands. Thus, regenera- 
tion method was not used as an independent variable. Sampling for the re- 
generation model was done plot by plot because of the heterogeneous nature 
of stands. Because the regeneration model predicts regeneration on a plot- 
by-plot basis, it is flexible for changes in forest management practices. The 
regeneration model should be applicable to a wide variety of silvicultural pre- 
scriptions, from stands that are treated homogeneously to stands that are 
treated to create heterogeneous stand conditions. 

It may be necessary and desirable to update the regeneration model in 
the future. Updates would incorporate changes in regeneration technology, 
account for long-term changes in climate, and represent shifts in species 
composition (such as increasing proportion of western white pine resistant 
to white pine blister rust [Cronartium ribicola]). The ability to update the 
regeneration model will depend on historical records being kept today. 
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APPENDIX A: HABITAT TYPE GROUPINGS 


This appendix documents results of analyses to combine similar habitat types into groups for the 
purpose of predicting regeneration establishment. An analysis of variance was performed by habitat 
type series (Pseudotsuga menziesii, Abies grandis, Thuja plicata, Tsuga heterophylla, and Abies lasio- 
carpa) using Scheffe’s S test to group similar habitat types, using a 0.05 significance level. Scheffe’s 
test can be used with an unequal number of observations per habitat type, and it is robust regarding 
nonnormality and heterogeneity of variance (Kirk 1982). 

The results of these analyses, by themselves, did not determine the final groupings. There is varia- 
tion within habitat type in years since disturbance, aspect, elevation, geographic location, and so on 
that is not accounted for in the analyses. The results were useful in developing preliminary group- 
ings that were reviewed by forest ecologists. Their suggestions were incorporated to arrive at the 
final groupings. 

In the following tables, TPSP means number of trees per stocked plot and SPSP means number of 
species per stocked plot. See table 2, footnote 1, for definitions of species abbreviations used in nam- 
ing habitat types. The stocking mean is the ratio of plots that were stocked. Habitat types under- 
lined with the same line were not significantly different at the 0.05 level. Habitat types with fewer 
than 25 plots sampled were omitted from the analyses, but were later placed in the group thought to 
be most similar. 


Pseudotsuga menziesii Series 


Stocking mean= 0.73 0.72 0.51 042 0.33 033 0.30 0.28 0:25 0.25 
PSME/ VACA LIBO VAGL AGSP SYOR SPBE CARU CAGE PHMA SYAL 


TPSP mean=3.88 3.68 3.16 3.18 282 275 (1254 228 192 1.71 
PSME/ LIBO VACA VAGL CARU CAGE SYAL AGSP PHMA SYOR SPBE 


SPSP mean = 1.63 1.42 11.27 1.26 1.24 TAG toe 5) Pal 2h ale 07, 
PSME/ VACA LIBO AGSP VAGL CARU SYOR SYAL PHMA CAGE SPBE 


Final groupings: 1 = PSME/VAGL, LIBO, VACA 
2 = PSME/CARU, CAGE, AGSP, FEID 
3 = PSME/PHMA, ACGL 
4 = PSME/SYAL, SPBE, SYOR, ARUV, Misc. 


‘Average of only 11 stocked plots. 


Abies grandis Series 


Stocking mean= 0.85 0.77 0.51 0.50 042 0.39 0.38 0.38 0.33 
ABGR/ CLUN LIBO CLUN! XETE ASCA PHMA VAGL SPBE ACGL 

-XETE 
TPSP mean = 5.36 4.94 436 3.26 291 265 262 2.34 2.23 


ABGR/ CLUN LIBO CLUN! XETE PHMA VAGL SPBE ASCA ACGL 
-XETE 


SPSP mean = 1.92 1.77 1.50 148 1.48 1.33 1.30 1.28 1.20 


ABGR/ CLUN LIBO XETE VAGL CLUN! PHMA ACGL SPBE ASCA 
-XETE 
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Final groupings: 5 = ABGR/LIBO, CLUN-XETE 
6 = ABGR/XETE, VAGL, COOC, VACA 
7 = ABGR/CLUN (except CLUN-XETE) 
8 = ABGR/SPBE, ACGL, PHMA, ASCA, SETR 


‘Except CLUN-XETE. 


Thuja plicata Series 
Stocking mean = 0.61 0.56 
THPL/ CLUN ASCA 
TPSP mean = 9.64 5.41 
THPL/ CLUN ASCA 
SPSP mean = 1.75 1.48 
THPL/ CLUN ASCA 


Final groupings: 9 = THPL/AIl 


Tsuga heterophylla Series 
Stocking mean = 0.72 0.59 
TSHE/ CLUN ASCA 
TPSP mean = 10.66 5.30 
TSHE/ CLUN ASCA 
SPSP mean= 2.05 1.85 
TSHE/ CLUN ASCA 


Final groupings: 10 = TSHE/All 


Abies lasiocarpa Series 
Stocking mean = 
0.74 0.68 068 066 066 065 061 056 054 047 043 043 0.34 0.30 0.25 


TSME ABLA ABLA ABLA ABLA ABLA ABLA ABLA ABLA ABLA TSME ABLA ABLA ABLA ABLA 
XETE VACA ALSI MEFELIBO VASC VAGL XETE CLUN STAM CLUN CAGE CARU CACA ACGL 


TPSP mean = 
(44 6:75 656 5:87) 25.085 5.0508 4:72 924: 71 4:60 4:33) paseo b4 ene Con ello leom 


TSME ABLA ABLA ABLA ABLA ABLA ABLA ABLA ABLA ABLA ABLA TSME ABLA ABLA ABLA 
XETE MEFE CLUN VACA STAM LIBO VASC ALSI CARU VAGL XETE CLUN CACA CAGE ACGL 


SPSP mean = 
86 1:84 1:70) 2:65) 26225058) ol 5 1.46 1.40 36) 1:28) 51245 2122) sic 


ABLA TSME ABLA ABLA TSME ABLA ABLA ABLA ABLA ABLA ABLA ABLA ABLA ABLA ABLA 
ALSI XETE LIBO CLUN CLUN MEFE VAGL XETE STAM VACA VASC CACA CARU ACGL CAGE 
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Final groupings: 11 = ABLA/VAGL, VASC, VACA 


12 = ABLA/XETE, LIBO 
13 = ABLA/CLUN, GATR 


14 = ABLA/CAGE, CARU, ACGL, SPBE 


15 = ABLA/MEFE, ALSI, and 


TSME/CLUN, XETE, MEFE, STAM 


16 = ABLA/CACA, STAM, LUHI 


lAverage of only 13 stocked plots. 
2Average of only 9 stocked plots. 


Percentage of Stocked Plots, Average Trees Per Stocked Plot, and Average 
Number of Species Per Stocked Plot by Habitat Type Group 


Group 


OBNHHhWNe 


Habitat types 


PSME/VAGL, LIBO, VACA 
PSME/CARU, CAGE, AGSP, FEID 
PSME/PHMA, ACGL 

PSME/SYAL, SPBE, SYOR, ARUV, Misc. 
ABGR/LIBO, CLUN-XETE 
ABGR/XETE, VAGL, COOC, VACA 
ABGR/CLUN (except CLUN-XETE) 
ABGR/SPBE, ACGL, PHMA, ASCA, SETR 
THPL/All 

TSHE/AI 

ABLA/VAGL, VASC, VACA 
ABLA/XETE, LIBO 

ABLA/CLUN, GATR 

ABLA/CAGE, CARU, ACGL, SPBE 
ABLA/MEFE, ALSI 

TSME/CLUN, XETE, MEFE, STAM 
ABLA/CACA, STAM, LUHI 
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Ratio of 
stocked 
plots 


Trees per 
stocked 
plot 


Species per 
stocked 
plot 


1.36 
1.21 
1.14 
1.13 
1.81 
1.49 
1.48 
1.28 
1.74 
2.05 
1.44 
1.56 
1.65 
1.23 
1.66 


1.36 


APPENDIX B: COEFFICIENTS FOR EQUATIONS 


Terms and abbreviations appearing in these tables may need some explanation. Variables are de- 
fined in table 3. Habitat type abbreviations are defined in table 2, footnote 1. Optimum aspect is 
calculated using the procedure developed by Stage (1976). The poorest aspects are 180 degrees from 
the optimum. Amplitude for probability of stocking equations uses SLO = 0.3, TIME = 10 years, and 
the probability of stocking centered on 0.5. Amplitude was not calculated when the optimum aspect 
was not significant. NS stands for not significant at the 0.05 probability level. In this appendix foot- 


notes appear in front of the items they refer to. That position should help avoid confusion with the 
superscripts for “squared” in BAA?, TPSP?, and ELEV’. 
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Table 9—Probability of stocking equations for version 2 of the Regeneration Establishment Model. The form of the equation is 


P= (1+e-8X))—1 
Thuja 
plicatal 
Pseudotsuga Abies Tsuga Abies 
menziesii grandis heterophylla lasiocarpa All 
Variable series series series series roads 
a iuneeemee (0), geeks ase ee) ae YONA PS eR PB ew Be eB) 
1 Intercept —0.82988 —-2.42977 —6.21669 —0.43035 —2.10728 
2 COS(ASP)*SLO'?*TIME'2 .07406 39253 .58797 .24625 0) 
3 SIN(ASP)*SLO'2*TIME'2 —.06721 .11738 .00742 -.01938 0 
4 SLO'**TIME'2 —.02119 —.31847 15254 —.09997 0 
5 TIME’ 0 55349 0 0 0 
6 SQREGT .21368 0 51481 .23914 .42067 
7 SQBWAF .12913 0 .44986 .22470 .18714 
8 BWB4 0 0 0 -—.11159 —.35584 
9 ELEV —.02706 .07568 .28859 0 (0) 
10 ELEV? 0 —.00083 —.00395 0 0 
11 BAA 0 .01382 .00795 0 01113 
12 BAA? 0 —.00008 —.00004 0 0 
13 In(BAA) 0 0 0 .13678 0 
14 COS(ASP)*SLO*BAA 0 -.01346 -—.01790 0 0 
15 SIN(ASP)*SLO*BAA 0 —.01348 —.00600 0 0 
16 ‘NONE 0 0 0 0 No data 
17 MECH —.18027 -.17355 —.18546 .06815 No data 
18 BURN -.48567 -.17355 —.34635 —.34280 No data 
19 PLANT .74846 0 0 0 0 
Habitat type groups 
20 PSME/VAGL 1.14975 No data No data No data 0 
21 PSME/CARU .07020 No data No data No data 0 
22 PSME/PHMA —.11583 No data No data No data 0 
23 PSME/SYAL 0 No data No data No data 0 
24 ABGR/LIBO No data 54137 No data No data 59764 
25 ABGR/XETE No data —.05775 No data No data .59764 
26 ABGR/CLUN No data 0 No data No data 59764 
27 ABGR/SPBE No data -.17732 No data No data 59764 
28 THPL/AIl No data No data 0 No data 1.35479 
29 TSHE/AIl No data No data .35600 No data 1.50736 
30 ABLA/VAGL No data No data No data .27860 1.16959 
31 ABLA/XETE No data No data No data .18176 1.16959 
32 ABLA/CLUN No data No data No data 0 1.16959 
33 ABLA/CAGE No data No data No data —.76407 1.16959 
34 ABLA/MEFE No data No data No data .40383 1.16959 
35 ABLA/CACA No data No data No data —.83775 1.16959 
National Forests 
36 Bitterroot 1.07708 0 No data 0 0 
37 Panhandle 0 0 0 0) 0 
38 Clearwater 0 -.47893 0 0 0 
39 Colville 0 0 0 0 0 
40 Deerlodge .73060 No data No data 0 0 
41 Flathead .73060 0 No data 0 0 
42 Gallatin .73060 No data No data 0 0 
43 Helena .73060 No data No data 0 0 
44 Kootenai No data 0 0 0 0 
45 Lolo 1.07708 0 0 0 0 
46 Nez Perce 0 -.41795 0 0 0 
47 Payette .28613 —1.08325 No data 0 0 
48 Boise .80554 —1.08325 No data 0 0 
Chi-square 22.05 33.38 7.60 26.88 27.16 
Error mean square 0.9923 1.0100 0.9973 1.0038 1.0307 
Percent stocking 30.90 49.54 65.37 55.97 51.24 
Total No. of plots 2,369 3,373 3,289 2,414 683 
Optimum aspect 318° We? ile 356° NS 
Amplitude 0.09 0.34 0.47 0.21 —_ 


‘Represented as part of the constant term to avoid a singular matrix. 
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Table 10—Weibull distribution functions for determining the number of trees 
per stocked plot 


TPSP = B[-In(1-X) ]"”°+ 1 
where 


TPSP = number of trees per stocked plot 

X =a uniform pseudo-random number in the interval [0,1] 

B = EXP[1.79862 + 0.64299*COS(ASP)*SLO 
— 0.34931*SIN(ASP)*SLO — 2.18751*SLO — 0.03815*ELEV 
+ 0.27987*SQREGT + 0.15874*SQBWAF 
+ 0.07241 if Abies grandis series 
+ 0.85008 if Thuja/Tsuga series 
+ 0.64994 if Abies lasiocarpa series] 

C = EXP[-— 0.33367 — 0.00751*ELEV + 0.07164*SQREGT 
+ 0.06127*SQBWAF] 


Table 11—Equations for determining the number of species on stocked plots. The form of the equation is P = (1+e~-@4i*?)—". 
The equation for one species is conditional on there being at least two trees on the plot, while other equations are 


conditional on the number of trees being equal to or greater than the number of species being predicted 


Variable 
i (X) 


1 Constant 
2 COS(ASP)*SLO 
3 SIN(ASP)*SLO 
4SLO 
5 REGT 
6 BWAF 
7 ELEV 
8 BAA 
9 PLANT 
10 TPSP 
11 TPSP? 
12 In(TPSP) 


Overstory climax 


13 'PSME 
14 ABGR 
15 THPL 
16 TSHE 
17 ABLA 


Chi-square 

Error mean square 
2Percent occurrence 
No. plots 

Optimum aspect 
Amplitude 


‘Represented as part of the constant term to avoid a singular matrix. 


1 species 2 species 3 species 4 species 5 species 
(B) (B) (B) (B) (B) 
1.39959 —0.44188 —2.05274 —6.24551 —2.04323 
0 57776 .07506 —1.27740 0 
0 .29407 —.47116 —.32990 0 
0 —.12877 —.84038 .42524 0 
—.02217 .01657 .00475 .00137 0 
—.03940 —.00025 .07635 .07188 0 
.02139 —.01104 —.01859 0 -.05291 
.00216 0 0 0 0 
—.39803 0 .28128 .61675 0 
0 —.01013 .05466 .05678 .02188 
0 0 —.00047 —.00024 0 
—1.01790 0 0 0 0 
0 0 0 0 0 
—.69564 .43695 .67734 2.30190 0 
—.77642 .42663 .90042 2.60261 0 
—1.22760 .36357 1.29021 3.08950 0 
—1.05898 72147 .90065 2.15632 0 
14.94 7.01 20.35 13.82 13.29 
0.9900 1.0028 0.9975 0.9223 0.9402 
40.34 39.24 19.79 6.87 2.35 
4,370 4,370 3,269 2,678 2,213 
NS 272 279° 194° NS 
— 0.10 0.07 0.20 — 


6 species 


(B) 
5.93804 


oooo0oO0CoO0 0 


.02007 


ooo0o°o 


Percent occurrence for dichotomous dependent variables is the percentage of plots on which the event occurred. For example, 1,763 of 
4,370 stocked plots had only one species; therefore, percent occurrence is 40.34. 
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Table 12—Probability of occurrence for advance regeneration on stocked plots. The form of the equation is P = (1+e-**))-" 


Ponderosa Western Lodgepole Engelmann Douglas- 


Variable pine 
U (X) (B) 
1 Constant -2.07555 
2COS(ASP)*SLO -2.53235 
3SIN(ASP)*SLO -.59257 
4SLO -.55116 
5 ‘NONE 0 
6 MECH —.56829 
7 BURN —.72614 
8 ROAD No data 
9 BAA .02275 
10 BAA? -.00024 
11 In(BAA) 0 
12 ELEV 0 
13 ELEV? 0 
14 TIME 0 
15 REGT —.12626 
16 BWAF -.08573 
17 BWB4 0 
18 OVER 1.05233 
19 BOTTOM 0 
Habitat type groups 
20 PSME/VAGL 0 
21 PSME/CARU 0 
22 PSME/PHMA 0 
23 PSME/SYAL 0 
24 ABGR/LIBO —2.59212 
25 ABGR/XETE -2.59212 
26 ABGR/CLUN -2.59212 
27 ABGR/SPBE —.70368 
28 THPL/All —2.59212 
29 TSHE/All No data 
30 ABLA/VAGL No data 
31 ABLA/XETE No data 
32 ABLA/CLUN No data 
33 ABLA/CAGE No data 
34 ABLA/MEFE No data 
35 ABLA/CACA No data 
National Forests 
36 Bitterroot 0 
37 Panhandle 0 
38 Clearwater 0 
39 Colville 0 
40 Deerlodge 0 
41 Flathead No data 
42 Gallatin No data 
43 Helena No data 
44 Kootenai 0 
45 Lolo 0 
46 Nez Perce 0 
47 Payette 1.53214 
48 Boise 1.53214 
Chi-square 16.39 


Error mean square 0.7653 
Percent occurrence 2.81 


No. plots 


Optimum aspect 


Amplitude 


3,456 
193° 
0.37 


larch 
(B) 


-5.00929 
.25604 
1.59651 
-.17195 
0 
0 
No data 
No data 
0 


-10395 


OnNOoOOoOOoOOCOOCOC0C Oo 


.49348 


oooooo+o000o 


— 
n 
o 
o 
rs 
© 


No data 


15.67 
1.0364 
1.55 
4,182 
81° 
0.24 


pine 


(B) 


—7.88764 
0 
0 
-1.71254 
0 
0 
0 
No data 
0 
0 
.29813 
.04579 


.47999 


OnNOOOOC Oo 


1.32705 
1.32705 
0 
0 
0 
1.93198 
1.93198 
0 
0) 
1.32705 
0 


oooooo0oo0°ooo9°;]o 


22.65 
1.2620 
2.76 
5,900 
NS 


spruce 


(B) 


—12.22361 
0 
0 
0 
0 
—.75266 
—.95539 
—.14421 


ooooo0o0°co°oo°c]oe 


oooo0o°o 


.94810 

.94810 
0 
1.14841 
1.14841 


7.71 
0.9844 
2.93 
5,494 
NS 


‘Represented as part of the constant term to avoid a singular matrix. 


fir 


(B) 


0.69112 
—.22493 
—.46164 
1.39088 
0 
—.91351 

—1.59609 

—1.26675 

.01519 
—.00008 
0 
—.01378 


ooooooo0o0o9°9co0o°ooe 


13.73 
0.9917 
13.34 
6,191 
244° 
0.08 
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White 
pine 
(B) 


—1.87330 
—1.58932 
-3.78659 
-3.97804 

0 
—.89870 
—.87392 
No data 
.02520 
-—.00031 

0 
.03855 

0 
—.22848 


No data 
No data 
No data 
No data 
0 
—1.22480 
0 
-1.22480 
—.42966 
—1.22480 
No data 
—1.22480 
—.42966 
No data 
—.42966 
No data 


No data 
0 
0 
-1.13793 
No data 
0 
No data 
No data 
0 
-1.13793 
No data 
No data 
No data 


11.75 
1.0493 
3.63 
4,046 
247° 
0.55 


Grand 
fir 


(B) 


—2.84267 
0 
0 
1.36587 
0 
—.99446 
—1.95615 
—1.36593 
.00725 
—.00003 
0 
.16438 
—.00208 
0 
—.11820 
—.10925 
0 
1.05382 


Western 
hemlock 


(B) 


—1.83783 
3.89786 
—.41924 
—.03178 
0 

—1.15450 

—3.18414 

—1.46549 


.19186 


o-o0oo0ooqo0o0co0o 0 


No data 
No data 
No data 
No data 
No data 
No data 
No data 
No data 
No data 
0 

No data 
No data 
No data 
No data 
No data 
No data 


No data 
0 
0 
0 
No data 
No data 
No data 
No data 
0 
0 
No data 
No data 
No data 


26.19 
0.9538 
16.60 
994 
354° 
0.53 


Subalpine 
fir 


(B) 


-—14.92353 
.23280 
—.84457 
.08100 

0 
—.72264 
—1.43964 
No data 
.00530 

0 

0 
.40971 
—.00330 

0 
—.06030 
—.05421 

0 
1.66590 

0 


No data 
No data 
No data 
No data 
1.09193 
0 

0 

No data 
0 

1.24962 
1.69269 
1.69269 
1.69269 
0 

2.37496 
2.37496 


ooooooo0°oo°coo0eo°c°o 


30.01 
0.9852 
8.62 
4,813 
285° 
0.13 


Western 
redcedar 


(B) 


—0.30393 
67711 
—1.01033 
1.76319 
0 
—.59370 
—1.38510 
—1.57013 
.00575 
0 
0 
—.03900 
0 
—.08185 
0 
0 
0 
1.40441 
1.11039 


No data 
No data 
No data 
No data 
No data 
No data 
No data 
No data 
0 
.65724 
No data 
No data 
No data 
No data 
No data 
No data 


No data 
—1.15548 

0 

0 

No data 

No data 

No data 

No data 

0 

0 

0 

No data 

No data 


20.87 
1.0074 
16.29 
2,326 
304° 
0.30 


Table 13—Probability of occurrence for subsequent regeneration. The form of the equation is P= (1+e-*)))" 


Ponderosa Western Lodgepole Engelmann Douglas- White Grand Western Subalpine Western 


Variable pine larch pine spruce fir pine fir hemlock fir redcedar 
i (X) (B) (B) (B) (B) (B) (B) (B) (B) (B) (B) 
1 Constant -3.92859 -6.63750 -0.96370 -17.66672 -2.62358 2.41585 -4.21749 -12.30531 —-7.73260 -3.62791 
2COS(ASP)*SLO -.71162 1.42144 —.68409 1.52230 57720 -.01901 .77652 2.56422 .87127 1.57737 
3 SIN(ASP)*SLO = -.21418 .77954 —.18537 1.00295 -.02488 -1.38626 .87180 —.43928 .43295 .73797 
4SLO —.65401 -—.49261 -3.19743 —-2.05762 .27942 -1.83669 -.45507 -1.40517 -1.75130 -1.79633 
5 ‘NONE 0 0 0 0 0 0 0 0 0 0 
6 MECH .67888 .98315 .64040 1.15332 .22750 .15480 -18407 .22827 .22322 .63192 
7 BURN .74077 1.12968 .38122 1.09008 .44157 .02302 .04756 .62188 .13249 1.06553 
8 ROAD .85366 1.51324 .79282 1.47729 .60009 55174 .56605 .65690 .93574 1.44289 
9 BAA 0 —.00973 —.03176 0 -.00587 -.00963 -.00453 0 -.01123 —.00399 
10 In(BAA) —.33190 0 0 -.24045 0 0 0 0 0 0 
11 ELEV —.03701 12605 (0) 0 —.02291 -.04806 .13693 45343 .05264 —.02599 
12 ELEV? 0 —.00170 0 0 0 0 -.00175 —.00583 0 0 
13 In(ELEV) 0 0 0 CIE) (0) 0 0 0 0 0 
14 TIME’ .46035 0 0 0 0 58166 0 74555 0 .92474 
15 SQREGT 0 0 0 55794 61203 0 55587 0 .54995 0 
16 SQBWAF 0 0 0 .71901 55493 0 .41238 0 .48266 0 
17 BWB4 0 0 0 0 0 0 —.24985 0 0 0 
18 OVER 1.07701 1.10264 1.58270 1.30915 23505 1.15252 .27597 .77804 .57990 .99736 
19 PLANT 1.43545 1.18734 .86778 1.79176 .83804 1.80834 .63057 0 .68549 0 
20 BOTTOM (0) 0 .72550 0 0 0 0 0 0 0 
21 UPPER .34694 0) 0 0 0 0 0 0 0 0 
22 RIDGE 34694 0 0 0 0 0 0 0 0 0 
Habitat type groups 
23 PSME/VAGL 0 .66031 .40172 No data 0 No data Nodata No data No data No data 
24 PSME/CARU .65576 0 .40172 No data 0 No data Nodata No data No data No data 
25 PSME/PHMA .65576 0 —.77148 No data 0 No data No data No data No data No data 
26 PSME/SYAL .65576 0 —.77148 No data 0 No data No data No data No data No data 
27 ABGR/LIBO —.39399 .66031 0 34244 0 0 0 No data 0 No data 
28 ABGR/XETE (0) .66031 .67108 34244 0 0 -.21632 Nodata 0 No data 
29 ABGR/CLUN —.39399 0 0 -.79350 0 0 0 No data 0 No data 
30 ABGR/SPBE .65576 0 —.77148 -.79350 0 0 0 No data No data No data 
31 THPL/AII No data .74374 —.54311 0 0 .65722 .22883 Nodata 0 0 
32 TSHE/AIl No data .89597 —.08682 21151 0 1.03756 0 0 1.28599 —.64033 
33 ABLA/VAGL No data .73817 1.73217 .09543 -1.10811 Nodata -2.44230 No data 1.03643 No data 
34 ABLA/XETE No data .73817 .67108 .09543 -.18290 0 -2.44230 Nodata 1.03643 Nodata 
35 ABLA/CLUN No data .73817 —.63157 .09543  -.18290 0 -.81250 Nodata 2.44933 Nodata 
36 ABLA/CAGE No data No data 1.73217 .09543 -1.10811 Nodata -2.44230 Nodata 1.03643 Nodata 
37 ABLA/MEFE No data .73817 .43158 -55866  -1.10811 0 —2.44230 No data 2.44933 No data 
38 ABLA/CACA No data No data —.63157 55866 Nodata Nodata Nodata No data 2.44933  Nodata 
National Forests 
39 Bitterroot 2.48159 0 0 —.52561 .61952 Nodata -1.80142 No data 0 No data 
40 Panhandle 0 0 -1.41402 1.10987 0 0 0 0 0 0 
41 Clearwater 0 0 -1.41402 .97309 0 0 0 0 0 0 
42 Colville 0 0 0 0 61952 -1.54687 0 0 0 0 
43 Deerlodge 0 1.87274 0 31152 .61952 Nodata Nodata No data 0 No data 
44 Flathead No data 1.87274 0 31152 61952 -1.54687 -1.80142 Nodata 0 No data 
45 Gallatin No data No data 0 —.52561 61952 Nodata Nodata No data 0 No data 
46 Helena No data No data 0 —.52561 .61952 Nodata Nodata No data 0 No data 
47 Kootenai 2.48159 1.87274 0 2.31698 0 0 -1.10737 0 81375 0 
48 Lolo 1.87642 1.87274 0 0 0 -1.54687  -1.10737 0 81375 0 
49 Nez Perce 1.87642 0 -1.41402 91425 0 Nodata 0 No data 0 0 
50 Payette 2.48159 0 —1.41402 -—.39580 -.66350 Nodata -1.10737 Nodata 0 No data 
51 Boise 2.48159 Nodata 0 .74887 -.66350 Nodata -1.10737 Nodata 0 No data 
Chi-square 26.51 17.49 53.94 34.66 38.80 24.18 19.81 9.63 21.83 12.26 
Error mean square 1.0310 0.9899 1.0079 0.9609 0.9620 1.1176 1.0041 0.9923 1.0842 1.0015 
Percent occurrence 13.21 8.23 10.18 8.96 20.22 7.81 35.61 26.16 8.06 15.69 
No. plots 3,664 5,712 6,250 5,494 6,191 4,326 5,293 994 5,125 2,326 
Optimum aspect 197° 29° 195° 33° 258° 269° 48° 350° 26° 25° 
Amplitude 0.11 0.24 0.11 0.27 0.09 0.21 0.17 0.37 0.14 0.26 


‘Represented as part of the constant term to avoid a singular matrix. 
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Table 14—Probability of occurrence for excess regeneration. The form of the equation is P = (1+e°@*?)—" 


Ponderosa Western Lodgepole Engelmann Douglas- 


Variable pine larch 
i (X) (B) (B) 
1 Constant —18.89119 -15.25330 
2COS(ASP)*SLO 0 1.58775 
3 SIN(ASP)*SLO 0 .44217 
4 SLO 0 —1.50273 
5 ‘NONE 0 0 
6 MECH 0 .75803 
7 BURN 0 .67282 
8 ROAD 0 1.92937 
9 BAA 0 0 
10 In(BAA) 0 0 
11 ELEV .66069 .43849 
12 ELEV? —.00690 -.00513 
13 TIME 0 0 
14 REGT 0 0 
15 BWAF 0 0 
16 BWB4 0 0 
17 OVER .91178 1.28724 
18 BOTTOM 0 0 
Habitat type groups 
19 PSME/VAGL 0 1.05732 
20 PSME/CARU 0 0 
21 PSME/PHMA 0 0 
22 PSME/SYAL 0 0 
23 ABGR/LIBO -.80025 1.05732 
24 ABGR/XETE -—.80025 1.05732 
25 ABGR/CLUN -1.63194 .43319 
26 ABGR/SPBE —.80025 0 
27 THPLAAIl —1.90005 .43319 
28 TSHE/All No data 1.20503 
29 ABLA/VAGL No data 1.55633 
30 ABLA/XETE No data 1.55633 
31 ABLA/CLUN No data 1.55633 
32 ABLA/CAGE No data No data 
33 ABLA/MEFE No data 1.55633 
34 ABLA/CACA No data No data 
National Forests 
35 Bitterroot 0 0 
36 Panhandle 0 0 
37 Clearwater 0 0 
38 Colville 0 0 
39 Deerlodge 0 2.39893 
40 Flathead No data 2.39893 
41 Gallatin No data No data 
42 Helena No data No data 
43 Kootenai 0 2.39893 
44 Lolo 0 2.39893 
45 Nez Perce 0 0 
46 Payette 0 0 
47 Boise 0 No data 
Chi-square 0.52 8.96 
Error mean square 0.9178 1.1135 
Percent occurrence 1.47 2.78 
No. plots 3,664 5,712 
Optimum aspect NS 16° 
Amplitude _ 0.24 


pine 
(B) 
—2.56629 
1.07703 
—.20678 
—1.75710 
0 
.05479 
-.61241 
—.47796 
0 
—.36016 
—.21089 
.00224 


ooooooo0oo0o0°c°eo9°c]e 


20.94 
0.9397 
4.45 
6,250 
349° 
0.16 


spruce 


() 


fir 


(B) 


-26.32721 -2.00802 


0 —0.13516 
0 —.39445 
-1.95903 .85319 
0 0 
99041 -.41157 
54089 -.49398 
1.52932 .03980 
0 0 
0 0 
74546 -.03835 
-.00650 0 
.07038 .06521 
0 0 
0 0 
0 0 
1.22240 .73419 
0 0 
No data .78603 
No data .78603 
No data 0 
No data 0 
0 0 
0 0 
.92664 0 
0 -—.60779 
1.85920 0 
2.33304 .49310 
61242 -.63641 
61242 52193 
1.85230 0 
0 -—.63641 
1.85230 0 
~ 1.85230 Nodata 
0 1.31556 
0 0 
0 0 
0 0 
0 1.31556 
0 0 
0 0 
0 1.31556 
0 0 
0 1.31556 
0 0 
0 0 
0 0 
9.00 23.69 
1.3310 0.9934 
2.60 10.32 
5,494 6,191 
NS 251° 
— 0.06 


"Represented as part of the constant term to avoid a singular matrix. 
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White 
pine 
(B) 
—2.66011 

—.21036 
—2.17665 
—2.61597 


.05466 


o-oo0o0oo9co 


No data 
No data 
No data 
No data 
-.67461 
-.67461 
0 

—.67461 
0 

0 

No data 
—.67461 
-.67461 
No data 
0 

No data 


No data 
0 
—.56300 
-1.86484 
No data 
—1.86484 
No data 
No data 
0 
—1.86484 
No data 
No data 
No data 


2.04 
0.9455 
3.33 
4,326 
264° 
0.32 


Grand 
fir 
(B) 

6.14484 
1.47745 
.26161 
-.47692 

0 

—.42492 
—.83001 
-.21219 
—.00489 
.13567 
.27891 
—.00366 
0 
.05824 
.05026 
0 
.43489 
0 


No data 
No data 
No data 
No data 
0 
-1.14939 
-.55579 
-1.14939 
0 
0 
-2.43159 
—2.43159 
—.88282 
-2.43159 
—2.43159 
No data 


Western 
hemlock 


(B) 


—9.31954 
4.31678 
—.90111 
—.13394 
0 
—.26230 
—.56556 
—.18974 
0 
0 

.39648 
—.00559 
.08965 


No data 
No data 
No data 
No data 
No data 
No data 
No data 
No data 
No data 
0 

No data 
No data 
No data 
No data 
No data 
No data 


No data 
0 
0 
0 
No data 
No data 
No data 
No data 
0 
0 
No data 
No data 
No data 


18.43 
0.9774 
24.95 
994 
358° 
0.58 


Subalpine 
fir 
(B) 
-7.40720 
1.03636 
.28255 
—1.52018 
0 
—.50965 
—.81329 
.21813 
0 
0 
.05698 


No data 
No data 
No data 
No data 
0 
0 
0 
No data 
0 
2.10127 
1.80532 
1.80532 
3.37135 
0 
3.37135 
2.10127 


1.20278 
0 

0 

0 

1.20278 
1.20278 
1.20278 
1.20278 
1.20278 
1.20278 


Western 
redcedar 


(B) 


=1291i79 
1.96111 
—.08096 
-.67317 
0 
—.06477 
-18186 
53623 
0 
0 
—.03679 
0 
.07178 


No data 
No data 
No data 
No data 
No data 
No data 
No data 
No data 
0 

-.43570 
No data 
No data 
No data 
No data 
No data 
No data 


No data 
0 
0 
0 
No data 
No data 
No data 
No data 


Table 15—Weibull distribution functions for determining delay to germination for advance best 
trees. The form of the equation is DELAY = B[—-In(1—X)]"’° + 1. Categories left 
blank use the corresponding coefficient for no budworm. “X” is a random number 
drawn from a uniformly distributed pseudo-random number generator. BWB4 is 
the number of budworm years in the 5 years before the last disturbance. BAA is 
the plot residual overstory area, reset to 1.0 if 0.0 


Species 


Ponderosa 
pine 
Western 
larch 
Lodgepole 
pine 
Engelmann 
spruce 
Douglas- 
fir 
White 
pine 
Grand 
fir 
Western 
hemlock 
Subalpine 
fir 
Western 
redcedar 


None 


0-25 
B= 8.98612 
C= 1.06847 


B= 9.76822 
C= 1.15258 


B= 7.35888 
C= .91230 


B =13.99027 
C= 1.05130 


B =13.12102 
C= 1.04325 


B= 6.69983 
C= 1.26253 


B =11.26918 
C= 1.12214 


B =13.60459 
C= 1.26706 


B =21.96234 
C= 1.10127 


B =17.77938 
C= 1.33722 


BWB4 (No. years) and BAA (ft?/acre) 


26+ 
10.31266 
1.47040 


27.10024 
1.31969 


32.9558 1 
1.23054 


21.77936 
1.41622 


22.18654 
1.21565 


13.43118 
1.27930 


18.24566 
1.08280 


19.60548 
1.28744 


32.17673 
1.31702 


24.24134 
1.54066 


1-3 


0-25 


8.38820 
1.33278 


16.85625 
1.16545 


11.39304 
1.13703 


20.20766 
1.66050 
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26+ 


18.02341 
1.50228 


18.60110 
1.18251 


19.99433 
1.35577 


30.05268 
1.83254 


10.58744 
1.47050 


18.83308 
1.30720 


16.15405 
1.33030 


25.51236 
1.87195 


16.55466 
1.20740 


17.66096 
1.93848 


13.93679 
1.31178 


20.55520 
1.24878 


Table 16—Equations predicting heights of advance regeneration. The form of the equation is HT = e Fix) 


Ponderosa Western Lodgepole Engelmann Douglas- 


Variable 
i (X) 


1 Constant 
2 In(AGE) 
3 COS(ASP)*SLO 
4 SIN(ASP)*SLO 
5 SLO 
6 ELEV 
7 ELEV? 
8 BAA 
9 In(TPSP) 
10 'NONE 
11 MECH 
12 BURN 
13 ROAD 
14 BWB4 
15 BWAF 
16 BOTTOM 
17 LOWER 
18 MIDSLOPE 
19 UPPER 
20 'RIDGE 


Habitat type groups 


21 2Dry PSME 
22 “Moist PSME 
23 ABGR 

24 THPL/TSHE 
25 'ABLA 


F-ratio 

R-square 

No. trees 
Optimum aspect 
“Amplitude 


pine 


(B) 


—1.52746 


86174 
0 
0 
0 
0 
0 
0 


-16174 


oo 00°00 


.07136 


ooo0o°o 


ooo0o°o 


40.70 
0.4806 
136 
NS 


larch pine 
(B) (B) 
-1.65811 -1.12979 
1.15051 91553 
0 .85450 
0 51484 
0 61714 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
57.53 42.85 
0.4979 0.4459 
60 218 
NS 31° 
— 1.82 


spruce fir 
(B) (B) 
—2.19638 -1.21096 
1.12147 1.07594 
0 0 
0 0 
0 0 
0 —.05148 
0 .00047 
-.00227 0O 
0 .14822 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 .00895 
0 .40658 
0 .21454 
0 .21097 
0 0 
No data 17475 
No data .02195 
0 .15324 
0 14169 
0 0 
109.32 92.65 
0.5430 0.5031 
187 estitat 
NS NS 


"Represented as part of the constant term to avoid a singular matrix. 
2Dry Pseudotsuga menziesii (habitat type group 1). 


3Moist Pseudotsuga menziesii (habitat type groups 2, 3, and 4). 


‘Height at the optimum aspect divided by height at the poorest aspect. 
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White 
pine 


(B) 


0.05586 
.84765 
—.79566 
.39278 
—.68673 
—.02835 
0 
—.00382 


oooooo0o0°o0o°co°o 


No data 
No data 


Grand 
fir 


(B) 


—2.26623 
1.02574 
—.06822 

-15558 
—.11461 


213.38 
0.4936 
2,200 
114° 
Usui 


Western 
hemlock 


(B) 


—0.70341 
.78900 


ooo0°c°o 


—.00322 
.07903 
0 
—.09283 
—1.16818 
—.35348 
0 


oooo0o°o 


No data 
No data 
No data 
0 

No data 


24.46 
0.3690 
258 
NS 


Subalpine 
fir 


(B) 


—1.88857 
.86013 


ooo 


12095 
0 
—.17341 
—.09718 

.14310 
—.04682 


oooo0o°o 


No data 
No data 
No data 
No data 
0 


82.48 
0.4716 
655 
NS 


Western 
redcedar 


(B) 


2.51134 
.13744 
—.08761 
.71709 
—.54406 
—.20298 
.00263 
0 
18618 


oooo0oo0o 0c 0co0o°c°o 


No data 
No data 
No data 
0 

No data 


40.35 
0.3013 
663 
97° 
1.54 


Table 17—Weibull distribution functions for determining delay to germination for subsequent 
best trees. The form of the equation is DELAY = B[-In(1—X)]"’ + 1. Categories left 
blank use the corresponding coefficient for the 0-2 year budworm category. “X” is 
a random number drawn from a uniformly distributed pseudo-random number gen- 
erator. BWAF is the number of years since the last disturbance with budworm 


Years since disturbance and BWAF 


TIME = 2-7 yr TIME = 8-12 yr TIME = 13+ yr 
Species 0-2 yr 3+ yr 0-2 yr 3+ yr 0-2 yr 3+ yr 
Ponderosa B= 3.81610 5.74622 9.36345 
pine C = 3.01975 2.09376 1.80925 
Western B = 5.23792 7.38005 10.42350 
larch C = 3.11598 3.04038 2.87196 
Lodgepole B=5.33757 6.78727 9.45827 
pine C= 4.16994 3.59937 2.39138 
Engelmann B= 5.36466 5.92251 7.44468 8.62548 9.69507 8.76074 
spruce C = 2.89777 3.31325 2.80504 5.21269 3.27745 2.30270 
Douglas- B= 4.34376 4.01218 6.55916 5.86172 9.16226 10.61297 
fir C = 2.33194 2.16561 2.63560 2.16010 2.21663 2.44953 
White B = 3.52946 7.62339 12.79801 
pine C= 1.71621 2.72466 3.98359 
Grand B= 4.17909 3.67409 5.88262 6.12256 8.49857 8.74195 
fir C = 2.47058 2.19139 2.25957 2.38522 2.00065 2.03146 
Western B= 4.33094 6.30802 8.63060 
hemlock C= 1.97408 2.35053 2.00997 
Subalpine B= 3.45725 4.11555 6.34975 6.19124 8.65545 8.82962 
fir C = 2.06804 1.98007 2.79933 2.84149 2.28687 2.46665 
Western B= 4.16284 7.52536 10.20937 
redcedar C= 2.03892 3.12279 2.68340 
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Table 18—Equations predicting heights of subsequent regeneration. The form of the equation is HT = e @:*) 


Ponderosa Western Lodgepole Engelmann Douglas- White Grand Western Subalpine Western 


Variable pine larch pine spruce fir pine fir hemlock fir redcedar 
i (X) (B) (B) (B) (B) (B) (B) (B) (B) (B) (B) 
1 Constant 0.07492  -1.47708 -0.66776 -2.79250 -2.21695 -1.86747 -2.51730 -2.52415 -2.29466  -1.40304 
2 In(AGE) 1.43793 1.19295 1.25739 1.40463 1.22826 1.27309 1.07561 1.41312 1.11917 .95548 
3 COS(ASP)*SLO 0 0 —.21656 0 -.15794 0 —.16952 0 0 —.56659 
4 SIN(ASP)*SLO 0 0 .03283 0 -.21290 0 .04787 0 0 .24004 
5 SLO 0 0 —.49383 0 -.09355 0 —.20352 0 0 .76533 
6 ELEV 0 0 —.03714 0 0 0 .03880 0 0 -.01515 
7 ELEV? 0 0 .00026 0 0 0 —.00055 0 0 0 
8 In(ELEV) —.61898 0 0 0 0 0 0 0 0 0 
9 BAA -.00549 -.00463 -.00848 -.00298 -.00370 -.00323  -.00421 -.00325 -.00434 —.00305 
10 In(TPSP) .07678 .08910 .13010 .07743 .13622 0 .19073 .15078 .16980 .16188 
11 ‘NONE 0 0 0 0 0 0 0 0 0 0 
12 MECH .20702 0 12941 -11501 .08694 0 -.09217 0 0 —.12496 
13 BURN 22443 0 .04401 .30555 .25036 0 .03063 0 0 .02469 
14 ROAD .14508 0 .12617 .08946 13741 0 -.12040 0 0 .15949 
15 BWAF 0 0 .01327 0 0 0 0 0 0 0 
16 BOTTOM 0 0 .36066 0 -.29274 0 -.07201 0 0 0 
17 LOWER 0 0 .20156 0 -.20692 0 —.15738 0 0 0 
18 MIDSLOPE 0 0 .30022 0 -.10117 O —.11644 0 0 0 
19 UPPER 0 0 .34617 0 -.24640 0 —.13008 0 0 0 
20 'RIDGE 0 0 0 0 0 0 0 0 0 0 
Habitat type groups 
21 ?Dry PSME 0 —.00159 —.07687 No data .01435 Nodata Nodata No data No data No data 
22 *Moist PSME 0 .01830 -.19234 Nodata -.19420 Nodata Nodata No data No data No data 
23 ABGR 0 .10002 .02146 22231 .21870 0 0 No data No data No data 
24 THPL/TSHE 0 .29961 —.05622 .26018 .40546 0 0 0 No data 0 
25 'ABLA 0 0 0 0 0 0 0 No data 0 No data 
F-ratio 149.59 78.08 57.24 128.43 92.98 159.00 212.47 180.03 127.53 41.94 
R-square 0.6116 0.4922 0.5719 0.6282 0.4943 0.4739 0.5021 0.6102 0.4272 0.4802 
No. trees 673 572 878 617 1,635 356 3,177 349 517 465 
Optimum aspect NS NS 171° NS 233° NS 164° NS NS 1572 
‘Amplitude — — 1.14 — 1.17 — euAl — — 1.45 


‘Represented as part of the constant term to avoid a singular matrix. 
2Dry Pseudotsuga menziesii (habitat type group 1). 

3Moist Pseudotsuga menziesii (habitat type groups 2, 3, and 4). 
‘Height at the optimum aspect divided by height at the poorest aspect. 
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Table 19—Weibull distribution functions for determining 
heights of excess trees. The form of the equation 
is HT = B[An(i-X)]"° +1. “X” is arandom num- 
ber drawn from a uniformly distributed pseudo- 
random number generator 


Years since last disturbance 


Species 2-7 yr 8-12 yr 13+ yr 
Ponderosa B = 2.77922 9.03331 14.13121 
pine C= .89932 1.07493 .93070 
Western B = 6.64373 11.42298 19.61887 
larch C= .90291 1.16616 1.30638 
Lodgepole B = 7.36042 10.92885 25.21441 
pine C = 1.14808 1.23233 1.11702 
Engelmann B= 1.46615 5.15927 9.27278 
spruce Cz= .72253 .73903 1.12551 
Douglas- B = 3.81608 8.16147 10.98770 
fir Cz= .99673 84541 .94804 
White B= 2.12146 5.06040 5.97955 
pine C= .74585 .78217 .84217 
Grand B = 3.08957 5.83018 10.10575 
fir C= .80068 83228 .95408 
Western B = 3.34771 6.80682 13.55346 
hemlock C= .56777 .89463 1.21404 
Subalpine B = 2.92136 4.58138 10.33328 
fir C= .88514 87156 1.04376 
Western B = 3.16951 4.50640 8.94054 
redcedar C= .64055 81354 .94349 


Table 20—Equations predicting the probability of site preparations. 
Only stands that were not deliberately site prepared 

were uSed in these analyses. The form of the equation 
is P= (1+e-@8/%))-" 


Variable Mechanical Burn None 
i (X) (B) (B) (B) 
1 Constant —-1.85203 -15.19530 1.04315 
2 COS(ASP)*SLO .49267 .05195 —.22095 
3 SIN(ASP)*SLO .19202 —.61358 .36958 
4 SLO —.96667 —.08902 .76911 
5 In(BAA) —.08592 —.37792 .26018 
6 ELEV .02494 53037 —.02969 
7 ELEV? 0 —.00491 0 
Habitat type series 
8 'PSME 0 0 0 
9 ABGR —.22684 .08709 .08573 
10 THPL —.60584 —.09800 .15116 
11 TSHE —.83269 —.75666 .68076 
12 ABLA —.30560 —.26382 .20339 
Chi-square 11.79 30.11 14.51 
Error mean square 0.9926 1.1297 1.0019 
Percent occurrence 18.66 9.40 66.00 
No. plots 6,024 6,024 6,024 
Optimum aspect eis 275° 121° 
Amplitude 0.08 0.09 0.06 


‘Represented as part of the constant term to avoid a singular matrix. 
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